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Abstract 

This  report  examines  the  structure  and  physical  properties  of  the  surface 
of  ice:  a  liquid-like  layer  on  the  ice  surface,  its  thickness  and  molecular 
structure,  surface  conductivity,  surface  electric  charge,  surface  potential, 
surface  diffusion  and  so  foith.  The  author  reviews  experimental  results  and 
theoretical  models  on  ice  adhesion,  gas  adsorption  and  ice  sintering.  Spe¬ 
cial  attention  is  given  to  the  results  of  studies  from  the  last  two  decades. 
Among  the  experimental  techniques  under  consideration  are  nuclear  mag 
netic  resonance.  X-ray  diffraction,  optical  ellipsometry,  measurements  of  ice 
surface  conductivity  and  dielectric  permittivity,  frictional  electrification  of 
snow  and  ice,  proton  channeling  and  others.  In  conclusion,  the  theoretical 
models  of  the  ice  surface  stri '  are  and  the  results  on  molecular  dynamics 
are  considered. 


For  conversion  of  SI  melric  units  to  U.S, /British  customary  units  of  measurement 
consult  Standard  Practice  for  Use  of  the  International  System  of  Units  (SI),  ASTM 
Standard  E380-89a,  published  by  the  American  Society  for  Testing  and  Mater¬ 
ials,  1916  Race  St.,  Philadelphia,  Pa.  19103. 

This  report  is  printed  on  paper  that  contains  a  minimum  of  50%  recycled 
material, 
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FOREWORD 


At  the  present  time,  thousands  and  thousands  of  people  around  the  world  deal  with  ice, 
snow  and  permafrost.  They  are  scientists,  educators,  engineers,  navigators,  meteorologists 
and  others.  While  a  small  fraction  of  these  people  contribute  to  the  knowledge  base  in  ice 
physics,  all  of  them  use  knowledge  from  it  frequently.  Moreover,  successful  applied 
research  is  based  upon  fundamental  science — one  more  reason  for  ice  specialists  to  have  a 
textbook  on  ice  physics  on  their  desks. 

The  first  modern  ice  physics  text  was  Fletcher's  book  on  The  Chemical  Physics  of  Ice  (1970). 
Fletcher's  book  is  in  typical  textbook  format:  it  is  reasonably  brief  and  easy  to  understand. 
He  touched  on  a  few  of  the  most  important  topics,  but  not  all  of  them. 

The  most  comprehensive  book  on  ice  physics  to  date  was  published  by  Hobbs  in  1974. 
Hobbs  considered  almost  all  of  the  basic  aspects  of  ice  as  understood  at  that  time.  Moreover, 
he  described  and  compared  se  /eral  (sometimes  opposing)  viewpoints.  This  fundamental 
and  rather  large  (837  pages)  book  is  commonly  known  as  the  "Ice  Bible"  by  specialists  in  the 
field.  In  1974  and  1975,  two  CRREL  Monographs  on  ice  W'ere  produced  by  John  Glen.  These 
were  briefly  arid  clearly  written  and  reviewed  almost  all  ice-related  subjects.  This  work  was 
(and  in  some  respects  still  is)  a  magnificent  introduction  to  ice. 

Finally,  in  1981  Maeno  wrote  a  simple,  popular  book  for  the  express  purpose  of  attracting 
people's  attention  to  the  subject. 

During  tire  past  20  years,  a  significant  amount  of  new  experimental  and  theoretical  work 
has  appeared,  dramatically  changing  our  views  on  ice  physics.  As  a  result,  we  are  now  able 
to  formulate  physical  laws  using  more  simple  and  direct  methods.  We  have  found  some  of 
the  physical  models  used  in  the  past  to  be  completely  wrong.  The  physics  of  ice  is  a  much 
better  developed  subject  than  it  was  20  years  ago. 

For  the  above  reasons,  I  feel  the  time  is  ripe  for  a  contemporary  book  on  ice  physics, 
incorporating  the  known  and  proven  with  almost  20  years  worth  of  material  not  covered 
by  previous  works. 

I  have  tried  to  prepare  a  "readable"  book,  and  not  one  that  requires  the  reader  to  be  a 
uniquely  educated  person.  It  is  my  intent  to  present  the  material  in  such  a  way  that  any 
reader  attracted  by  the  title  Ice  Physics  will  be  able  to  comprehend  it.  This  is  quite  difficult 
for  a  book  dedicated,  not  to  a  particular  field  of  knowledge,  but  to  a  specific  material. 
Indeed,  for  ice  it  means  we  have  to  consider  a  wide  variety  of  subjects,  including  quantum 
chemistry,  solid  state  physics,  tire  theory  of  elasticity,  ionic  conduct!  vity,  synchrotron  X-ray 
topography,  crystal  growth,  the  physics  of  surfaces  and  more 

The  primary  goal  is  to  produce  as  simple  a  book  as  possible  without  sacrificing  scientific 
accuracy.  Experimental  facts,  physical  ideas  and  theories  will  be  strongly  organized  and 
bound  together  cohesively.  Tire  reader  will  bi  introduced  to  a  wide  variety  of  material  on 
a  step  by  step  basis.  Then  the  picture  will  be  /hole. 

To  accelerate  publication,  this  book  will  appear  first  in  tire  form  of  a  series  of  joint 
CRREL-Dartmouth  reports,  later  to  be  published  in  CRREL/s  Monograph  series,  on: 

1 .  The  structure  of  ordinary  ice 

Part  I:  "Ideal"  structure  of  ice.  Ice  crystal  lattice 
Part  II:  Defects  in  Ice 
Volume  1:  Point  defects 
Volume  2:  Dislocations  and  planar  defects 

2.  Electrical  properties  of  ice 

Part  I:  Conductivity  and  dielectric  permittivity  of  ice 
Part  II:  Advanced  topics  and  new  physical  phenomena 


3.  Optical  properties 

4.  Electro-optical  effects  in  ice 

5.  Thermal  properties 

6.  Mechanical  properties  of  ice.  Elasticity  and  anelastic  relaxation.  Plastic  proper¬ 

ties.  Fracture  of  ice 

7.  Electromechanical  effects  in  ice 

8.  Hie  Surface  of  ice 

Part  I:  Experimental  results  on  the  structure  and  properties  of  the  ice  surface 
Part  II:  Theoretical  models  of  the  ice  surface 

9.  Other  forms  of  ice  and  their  properties 

10.  Ice  in  space 

11.  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  convenient  to  the  author.  The  present 
report  is  the  fifth  in  the  series. 


NOMENCLATURE 

a  length  of  hop 

mTi  length  of  screening  by  Bjerrum  defects 

az~2  length  of  screening  by  ions 

aT1  screening  length  when  I  e,<ps  I  »  kaT 
OLo  fraction  of  surface  inoleailes  in  the  liquid¬ 

like  layer  with  dipole  moments  (i.e.,  pro¬ 
tons)  oriented  towards  the  surface 
B  magnetic  induction 

Cs  concentration  of  adsorbed  molecules  on 

the  surface 

Cm  maximum  concentration  of  adsorbed 

molecules 

d  thickness  of  the  quasi-Uquid  layer 

Di  diffusion  coefficient 

A,v|/,d  angles 
Ea  heat  of  adsorption 

E,  E  electric  field  strength 

e;  electric  charge  of  i'h  type  of  charge  carrier, 

i  =  1,2,3, 4 

Eo  dielectric  permittivity  of  vacuum 

Es  static  dielectric  permittivity  of  ice,  c$  ®  100 

(T  =  -10°C) 

Es  activation  energy  of  surface  conductivity 

E '  effective  electric  field  acting  on  a  molecule 

e„  high-frequency  dielectric  permittivity  of 

ice,  Eu*  =  3.2 
d>  =  3.85  ksTrM. 

$  wetting  angle 

G  Gibbs  free  energy 

g  gyromagnetic  constant 

y  surface  energy 

Ysl  surface  energy  of  ice/  water  interface 

Yvl  surface  energy  of  vapor/ water  interface 

Yvs  surface  energy  of  vapor/ ice  interface 

I  current 

t|  viscosity 


tp  electric  potential 

<ps  surface  electric  potential 

k  absorption  coefficient 

ko  Boltzman  constant 

L  thickness  of  the  surface  layer 

"ksXsa  surface  charge  density 

Xsa  surface  charge  absorbed  on  the  ice  inter¬ 
face 

m  dipole  momentum  of  a  water  molecule 
«  refractive  index 

ne  concentration  of  Bjerrum  defects 

«H 20  water  molecule  concentration  in  ice 

«i  charge  carrier  concentration,  i  =  1,2,3, 4 

ti\  refractive  coefficient  of  ice 

hw  refractive  coefficient  of  water 

Wi0„  concentration  of  ions 

V  frequency  of  ion's  hops 

P  pressure 

P  polarization  vector 

Pc  electronic  polarization  vector 

Pm  molecular  polarization  vector 

Q  electric  charge 

P  'll  Al 

p(x)  density  of  the  space  charge 

Pl  water  density 

r±  reflection  coefficient  for  light  polarized 

perpendicular  to  the  plane  of  incidence 
rj|  reflection  coefficient  for  light  polarized 

parallel  io  the  plane  of  incidence 
ru  constant 

s  surface  area 

Ob  Bjerrum  defect  conductivity 

Ob  bulk  conductivity 

<7jon  ionic  conductivity 

cf0  low-frequency  bulk  conductivity 

0S  surface  conductivity,  f2-1 


specific  bulk  conductivity  of  the  ice  sur  ¬ 

^MS 

Maxwell  relaxation  time 

face  layer,  il'hn1 

T 

1  m 

melting  point  of  ice 

adhesive  strength 

V 

velocity 

entropy  of  fusion  per  unit  mass 

X 

coordinate 

high-frequency  bulk  conductivity 

CO 

angular  frequency 

temperature 

co0 

resonance  frequency 

time 

Wi 

electron  work  function  of  ice 

correlation  time  of  molecular  motion 

wm 

electron  work  function  of  a  metal 

mean  lifetime  of  a  molecule  adsorbed  on 
the  surface 

M 

chemical  potential 

Debye  relaxation  time 

Debye  relaxation  time  of  the  surface  layer 

Mi 

chargecarriermobility,  /=l,2/3,4forH30+, 
OH-  ions,  D-  and  L-defects  respectively 

of  ice 
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VICTOR  F.  PETRENKO 

INTRODUCTION 


Since  childhood,  we  have  accustomed  our¬ 
selves  to  matter  being  either  solid  or  liquid  or  gas. 
Later,  we  become  acquainted  with  the  fourth  kind 
of  matter  called  plasma,  which  is  very  rarefied  and 
has  an  extremely  large  electrical  conductivity. 

Can  you  imagine  a  simultaneous  combination 
of  certain  properties  of  solid,  liquid  and  plasma 
matter?  It  is  present  on  the  surface  of  the  ice  upon 
which  we  walk  each  winter,  unsuspecting  that  it 
possesses  such  unique  characteristics. 

When  we  bump  into  a  sharp  comer  or  hurt 
ourselves  with  a  needle  or  cut  a  finger  with  a  knife, 
we  confirm  that  surfaces  of  solids  are  indeed  solid. 
This  is  true  in  so  many  cases  that  it  necessarily 
becomes  a  rule. 

In  the  case  of  ice,  however,  some  simple  obser¬ 
vations  make  us  doubt  that  it  has  a  rigid  and  per¬ 
manent  surface.  For  example,  if  we  bring  two  ice 
balls  together,  then  right  before  our  eyes,  a  small 
"neck"  starts  to  grow  between  them  (Fig.  1).  In 
time,  further  growth  will  lead  the  two  balls  to 
become  one.  In  this  simple  experiment,  some  sub¬ 
stance  can  be  considered  to  flow  between  the  balls 
thiough  the  "neck."  Since  ice  in  bulk  form  is  ap¬ 
parently  solid  (have  you  ever  fallen  on  ice?),  it  is 
most  likely  the  surface  that  flows. 


a.  On  contact.  b.  319  minutes  later. 

Figutel.  Growth  of  a  neck  between  two  ice  spheres  0.5 
mm  in  diameter  at-6"C  (from  Hobbs  1 974). 


In  a  process  similar  to  the  ice  balls'  unification, 
airy,  fluffy  snowflakes,,  after  lying  together  for 
some  time,  sinter  first  into  dense,  heavy  snow,  and 
next  into  solid  ice.  It  is  this  process,  going  on  for 
thousands  and  thousands  of  years,  that  has  formed 
the  kilometer-thick  ice  sheets  in  Antarctica  and 
Greenland  (mass  transport  through  the  vapor 
phase  and  a  plastic  flow  are  also  involved  in  this 
sintering  process). 

Evidently,  the  first  documented  guess  for  the 
existence  of  a  special  layer  on  the  surface  of  ice  was 
the  work  by  Faraday  (1859).  On  the  grounds  of  his 
long-standing  experiments  on  the  adhesion  of  ice, 
Faraday  called  this  layer  a  "liquid-like"  layer.  An¬ 
other  common  expression  is  "quasi-liquid"  layer. 
Now,  the  existence  of  some  special  layer  on  the  ice 
surface  has  been  proven  by  many  experimental 
techniques.  Hundreds  of  papers  are  dedicated  to 
the  investigation  of  the  structure  of  this  layer  and 
its  physical  properties.  Hie  most  essential  contri¬ 
butions  to  the  store  of  our  knowledge  about  this 
liquid-like  layer  were  made  during  the  past  two 
decades. 

In  the  present  report,  I  shall  describe  the  results 
of  the  experimental  research  that  provides  direct 
information  on  the  structure  and  thickness  of  this 
layer:  X-ray  diffraction,  proton  channeling,  ellips- 
ometry  and  NMR.  Then,  numerous  results  of  stud¬ 
ies  of  the  physical  properties  of  the  ice  surface  will 
be  considered,  An  analysis  of  these  properties  by 
means  of  comparing  various  models  also  allows  us 
to  obtain  indirect  information  about  the  structure 
and  thickness  of  the  special  layer  on  the  ice  surface. 
We  shall  not  consider  ice  friction  here,  since  sever¬ 
al  good  reviews  on  this  topic  have  been  published 
recently  (see,  for  instance,  Aharon  et  al.  [1988]  and 
Colbeck  [1992]).  Nevertheless,  when  all  the  reports 
are  united  into  a  book,  a  chapter  on  ice  friction  will 
be  included. 

In  conclusion,  major  known  theoretical  models, 
describing  the  structure  and  properties  of  the  liq¬ 
uid-like  layer,  will  be  discussed. 


PART  I.  EXPERIMENTAL  RESULTS  ON  THE  STRUCTURE 
AND  PROPERTIES  OF  THE  ICE  SURFACE 


X-RAY  DIFFRACTION 

When  considering  tire  structure  of  the  quasi¬ 
liquid  film  on  the  ice  surface,  we  must  first  ans¬ 
wer  the  question:  is  the  structure  of  the  film  or¬ 
dered  or  not?  That  is,  does  the  structure  have  a 
long  range  order?  An  i  if  so,  then  is  this  structure- 
different  from  that  of  ice?  It  seems  that  the  easiest 
way  to  answer  this  question  is  to  use  any  known 
X-ray  diffraction  technique.  However,  owing  to 
the  extremely  small  thickne- c  'f  the  film,  it  is  dif¬ 
ficult  to  differentiate  its  contmution  to  the  dif¬ 
fraction  from  that  of  the  ice  bulk.  I  am  aware  of 
only  one  such  attempt,  wnere  Kouchi  et  al.  (1987) 
investigated  X-ray  diffraction  at  the  surfaces  of 
polycrystalline  and  monocrystalline  ice  in  a  tem¬ 
perature  range  of  -0.5  to  -10°C.  Glancing  beams 
were  used  (angle  of  incidence  =  2°)  to  make  the 
contribution  of  the  surface  as  large  as  possible. 
The  Ka  line  of  the  Cu  X-ray  spectrum  was  used  in 
the  experiments. 

F:  '"'.re  2  snows  X-ray  diffraction  patterns  of 
tbit  surface  of  the  polycrystalline  ice  sample  ob¬ 
tained  at  -9.3,  -2.1  and  -1.1  °C.  At  T  =  -9.3  and 
-2.1°C,  the  diffraction  peaks  that  are  characteris¬ 
tic  for  the  crystal  are  clearly  visible,  although 
their  magnitude  decreases  at  -2.1°C.  At  T  -  -1.1 
°C,  those  peaks  virtually  vanish,  being  re- 


50  40  30  20 
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Figure  2.  Intensity  of  the  beam  diffracted  from  the 
polycrystalline  ice  surface  as  a  function  of  the  angle 
between  incident  and  reflected  beams  (after  Kouchi  etal, 
1987). 


placed  by  a  blur;_d,  circular  "halo"  cen¬ 
tered  at  2d  =  26°.  IT'S  is  a  clear  eviden  e 
that,  at  1 ...  cem,  e,  the  long- range 
order  in  the  arrangement  of  mole¬ 
cules  in  the  surface  layer  \ 

Since  the  intensity  of  t1'"  ’  *ic  00  £r* 

used  to  estimate  the  disorc  ’  -  " 

thickness,  it  is  interesting  to  to,.  me 
temperature  dependence  of  this  intensity, 
which  is  illustrated  in  Figure  3.  First,  we 
can  conclude  from  this  figure  that  the  lay¬ 
er's  thickness  is  independent  of  the 
crystallographic  orientation  of  the  sur- 


solid  :  polycrystal 
open  :  single  crystal 


\  V  (oooi) 


face.  Second,  the  layer  ceases  to  be  detect¬ 


able  by  this  method  at  T  <  -2°C.  Unfortu¬ 
nately,  Kouchi  e.  al.  did  not  attempt  to 


temperature  (°C) 


estimate  an  absolute  magnitude  of  the  Figure  3.  Changes  in  intensity  of  diffraction  halo  at  2d  =  28°  with 


layer's  thickness. 


temperature  (after  Kouchi  et  al.  1987). 
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PROTON  CHANNELING 

Another  interesting  experimental  method  of 
determining  the  thickness  and  structure  of  the 
quasi-liquid  layer  was  used  by  Golecki  and  Jac- 
card  (1977, 1978).  This  is  the  so-called  channeling 
effect,  discovered  by  Oen  and  Robinson  (1962). 
The  channeling  effect  manifests  itself  in  the  fact 
that,  if  a  beam  of  heavy  charged  particles  (usually 
protons  or  a  particles)  travels  parallel  to  some 
crystallographic  direction,  then  the  particles  of 
that  beam  experiencing  low-angle  scattering  pen¬ 
etrate  the  crystal  at  very  long  distances  and  can 
stay  inside  the  channels  formed  by  atomic  rows. 
The  magnitude  of  this  effect  is  usually  expressed 
in  terms  of  the  ratio  of  the  flux  of  backscattered 
particles  when  the  direction  of  the  flux  coincides 
with  one  of  the  crystallographic  directions  to  that 
when  the  direction  is  arbitrary. 

The  backscattering  can  take  place  only  at  nu¬ 
clei  of  atoms  that  are  much  heavier  than  the  parti¬ 
cles  in  the  beam  and  resembles  the  classic  Ruther¬ 
ford  scattering.  The  diffraction  of  such  particles  at 
the  crystal  lattice  is  insubstantial,  since  the 
characteristic  de  Broglie  wavelength  is  several 
orders  of  magnitude  less  than  the  lattice  constant. 

Most  of  the  backscattered  particles  experience 
only  a  single  strong  interaction  with  a  lattice 
atom.  Using  the  theory  of  Rutherford's  scatter¬ 
ing,  we  can  calculate  the  energy  of  the  reflected 
particle  as  a  function  of  the  scattering  angle.  De¬ 
partures  from  this  predicted  value  contain  infor¬ 
mation  about  the  depth  at  which  the  scattering 
took  place  (more  exactly,  about  the  path  the  parti¬ 
cle  has  traveled  back  and  forth  within  the  matter). 
The  loss  cf  the  kinetic  energy  of  the  particles  is 
attributable  to  their  interaction  with  tire  valence 
electrons  and  is  almost  independent  of  the  direc¬ 
tion  of  the  motion.  Having  "calibrated"  such  loss¬ 
es  for  a  particular  material  (or  having  calculated 
them  theoretically),  we  can  establish  the  depth  at 
which  the  scattering  took  place  with  a  rather  high 
accuracy. 

The  channeling  is  extremely  sensitive  to  the 
degree  of  perfection  of  a  crystal  structure  and 
therefore  can  be  efficiently  used  for  measure¬ 
ments  of  the  magnitude  of  thermal  vibrations,  the 
degree  of  disordering,  the  concentration  of  point 
defects  and  so  forth.  One  of  the  comprehensive 
reviews  on  the  subject  is  that  by  Gemmel  (1974). 

In  1977  and  1978  this  effect  was  explored  by 
Golecki  and  Jarcard  for  the  investigation  of  the 
(0001)  surface  of  ice  crystals.  The  authors  mea¬ 
sured  the  backscattering  (angle  150°)  of  100-keV 
protons  from  the  surface  of  monocrystalline  ice  in 


a  wide  temperature  range  from  -2  to  -130°C. 
They  found  that,  at  temperatures  higher  than 
-40°C,  the  surface  layer  is  characterized  by  much 
higher  disordering  in  relation  to  the  bulk.  Analy¬ 
sis  of  the  behavior  of  the  disordering  process 
prompted  Golecki  and  Jaceard  to  conclude  that 
the  disordering  is  caused  by  a  significantly  high¬ 
er  magnitude  of  thermal  vibrations  of  oxygen 
atoms  near  the  surface,  rather  than  the  absence  of 
symmetry  in  the  structure.  Thus,  at  T  =  -1.8°C  the 
average  amplitude  of  the  thermal  vibrations  of 
oxygen  atoms  near  the  surface  exceeds  that  of  the 
bulk  by  3.3  times.  Because  of  that,  Golecki  and 
Jaccard  suggest  not  using  the  term  'liquid-like' 
layer  at  all. 

To  estimate  the  thickness  of  the  special  surface 
layer  of  ice,  Golecki  and  Jaccard  calculated  an 
equivalent  thickness  of  an  amorphous  layer  (to¬ 
tally  disordered)  that  would  result  in  the  same 
magnitude  of  scattering.  The  results  of  such  cal¬ 
culations  are  depicted  in  Figure  4.  It  is  worth 
noting  chat  the  half-thickness  of  the  layer  t/2  ob¬ 
tained  is  extremely  large.  The  temperature  de¬ 
pendence  of  the  thickness  t  was  found  to  be 

d(nm)  =  (94  ±  17)  -  (54  ±  14)  log(273  -  T).  (1) 


273  -T(K) 

Figure  4.  Equivalent  "amorphous"  thickness  of  the 
disordered  region  at  the  ice  surface,  deduced  from  the 
measurements  with  the  aid  of  Meyer's  (1971)  theory. 
The  solid  line  corresponds  to  least-squares  fit.  The 
density  of  the  surface  region  was  assumed  equal  to  that  of 
ice  (0.917  g  enr3)  (after  Golecki  and  jaccard  1 978). 
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OPTICAL  ELLIPSOMETRY 

Ellipsometry  is  an  optical  reflectance  tech¬ 
nique  in  which  the  change  of  the  state  of  polariza¬ 
tion  of  light  upon  reflection  from  a  surface  is 
measured.  Even  though  the  method  dates  back  to 
Drude  (1889),  ellipsometry  has  recently  attracted 
considerable  attention  because  of  a  wide  variety 
of  applications  >n  modem  solid  state  optics,  espe¬ 
cially  in  the  physics  of  surfaces,  interfaces  and 
thin  films.  In  the  study  of  the  ice  surface,  optical 
ellipsometry  was  used  first  by  Beaglehole  and 
Nason  (1980),  and  later  by  Furukawa  et  al. 
(1987a,b). 

The  change  of  the  state  of  polarization  upon 
reflection  can  be  expressed  in  terms  of  the  ratio  of 
two  complex  reflection  coefficients  and  r±  for 
light  polarized  parallel  and  perpendicular  to  the 
plane  of  incidence,  r  being  the  ratio  of  reflected 
and  incident  electric  field  strength.  The  complex 
quantity 

p  =  tany  exp  (z'A)  (2) 

defines  the  two  ellipsometric  angles  A  and  vy.  The 
two  angles  completely  determine  the  two  optical 
constants  n  (refractive  index)  and  k  (absorption 
coefficient)  of  an  isotropic  reflecting  medium. 
Mathematical  dependence  of  n  and  k  on  n  and  vp, 
calculated  from  Fresnel's  formula,  was  found,  for 
example,  by  Bootsma  and  Meyer  (1969). 

The  application  of  ellipsometry  in  surface 
physics  is  favored  because  in  suitable  cases  a  lay¬ 
er  with  average  thickness  as  small  as  Vioo  of  a 
monolayer  can  be  detected,  if  in  the  layer  n  and  k 
are  different  from  that  of  the  bulk. 

Ellipsometnc  data  are  usually  analyzed  with 
computers  because  the  two  measured  variations 
8A  and  5 \\i  are  not  sufficient  to  determine  all  pa¬ 
rameters  (optical  constants,  thickness)  even  for  a 
single  surface  film. 

Beaglehole  and  Nason  (1980)  employed  ice 
monocrystals  grown  following  the  Bridgeman 
method  using  twice-distilled  water.  The  ice  sur¬ 
face  with  the  orientation  (0001)  (basal  plane)  or 
(1010)  (prismatic  plane)  was  milled  with  a  warm 
razor  immediately  before  ellipsometrical  tests. 
Water  forming  simultaneously  on  the  ice  surface 
was  removed  by  the  same  razor,  resulting  in  a 
surface  of  unruffled  cleanness. 

Measurements  in  the  0  to  -30'JC  temperature 
range  revealed  a  considerable  distinction  in  the 
behavior  of  the  surface  layer  at  the  basal  and  pris¬ 
matic  planes.  Assuming  that  the  refraction  coeffi¬ 
cient  of  the  surface  layer  is  the  same  as  that  of 


water,  Beaglehole  and  Nason  found  that  the  film's 
thickness  on  the  prismatic  surface  ( d )  was  about 
130  A  at  T  =  -1.5°C,  decreasing  slowly  to  a  mono- 
atomic  thickness  at  T  =  -10°C.  They  found  a  sur¬ 
face  layer  on  the  basal  plane  only  at  temperatures 
greater  than  -0.3°C. 

Furakawa  et  al.  (1987a,b)  significantly  ad¬ 
vanced  use  of  the  ellipsometry  technique  for  de¬ 
termining  the  transition  layer  thickness  on  the  sur¬ 
face  of  ice.  First,  they  used  improved,  more  perfect 
ice  crystals,  grown  from  twice-distilled  and  de¬ 
gassed  water.  The  degree  of  purity  of  ice  plays  an 
important  role,  since  even  a  small  salt  concen¬ 
tration  can  sizably  alter  the  thickness  of  the  transi¬ 
tion  (quasi-liquid)  layer  on  the  ice  surface.  Second, 
they  examined  atomically  flat  (0001)  and  (1010) 
surfaces  of  the  natural  facets  of  a  so-called  nega¬ 
tive  crystal.  A  negative  crystal  is  formed  by  pump¬ 
ing  water  vapor  out  of  a  small  void  inside  a  large 
monocrystal.  Besides,  they  made  their  ellipsomet¬ 
ric  measu.  ments  under  conditions  of  equilibrium 
between  ice  and  its  vapor. 

Finally,  Furakawa  et  al.  interpreted  their  data, 
not  using  the  assumption  that  the  refraction  coeffi¬ 
cient  of  the  surface  film  (n,)  is  the  same  as  that  of 
water  («w),  but  calculating  it  as  a  fitting  parameter. 
They  found  that  a,  =  1.330  for  both  faces  (0001)  and 
(101C),  which  is  Noser  to  water,  nw  =  1.3327,  than 
to  ice,  «i  =  1.3079.  This  implies  that  the  density  of 
tlr  layer  is  quite  close  to  that  of  water. 

Figure  5  shows  temperature  dependence  of  the 
transition  layer  thickness  as  found  by  Furakawa  et 
al.  from  their  ellipsometrical  measurements.  As 
follows  from  this  plot,  the  thickness  of  the  transi¬ 
tion  layer  drastically  rises  as  temperature  increas- 

Tempsrature  (°C) 


basal  (0001)  face 
prism.  (1010)  face 


Figure  5.  Temperature  dependence 
of  the  relative  phase  change  A  and 
the  thickness  of  liquid-like  layer  d 
(after  Furukawa  et  al,  1 987a), 
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es  above  -1°C,  exceeding  150  nm.  A  difference 
between  the  film's  thicknesses  on  the  (0001)  face 
and  the  (1010)  face  is  especially  noticeable  in  the 
temperature  range  from  -4  to  -1°C.  Optical  ellip- 
sometry  can't  detect  the  presence  of  a  liquid  film 
on  an  ice  surface  at  temperatures  below  -4°C. 

Very  recently,  Elbaum  et  al.  (1993)  performed 
optical  measurements  of  a  liquid-like  layer  using 
a  simplified  ellipsometry  technique,  which  was 
termed  Brewster  reflectometry.  The  idea  was  to 
measure  the  refractive  index  profile  near  a  plane 
interfac,"'  between  uniformmedia.  If  a  collimated 
Ugh:  beam  with  a  vector  E  lying  in  the  plane  of 
incidence  fails  at  the  Brewster's  angle,  then  tire 
intensity  of  tire  reflected  light  will  be  zero  for  an 
interface  between  two  homogeneous  media.  But, 
since  the  interface  has  a  finite  thickness,  the  re¬ 
flectivity  is  slightly  different  from  zero.  Elbaum 
et  al.  assumed  that  the  interface  is  composed  of 
water  vapor  or  its  mixture  with  air  having  the 
refraction  coefficient  n  =  1,  a  liquid  film  of  thick¬ 
ness  L  and  n  =1.333  (same  as  water)  and,  finally, 
ice  with  n  -  1.309.  Given  such  conditions,  tire 
problem  of  the  intensity  of  the  reflected  ligl  t  can 
be  solved  exactly,  and  they  used  this  solution  to 
calculate  the  thickness  of  the  liquid  film  on  the  ice 
surface  from  the  measurements  of  intensity  of  a 
polarized  light  reflected  at  the  Brewster's  angle. 
Though  this  method  has  a  lower  sensitivity  com¬ 
pared  to  con  -entional  ellipsometry  (L  2  12  A, 
Elbaum  et  al.  [1993]),  the  described  work  was  the 
first  where  the  effect  of  air  on  tire  thickness  and 
conditions  of  existence  of  the  quasi-liquid  layer 


T  (°C) 

Figure  6,  Collected  results  on  the  liquid-like  film  thick¬ 
ness  from  four  reflection  experiments.  The  inset  shows 
the  same  results  in  smaller  thickness  range  (after  Elbaum 
et  al.  1993). 


was  observed  and  investigated.  In  the  absence  of 
air,  the  quasi-liquid  film  could  not  be  detected  at 
the  prismatic  face  of  ice  within  the  experimental 
error  (12  A).  At  the  same  time,  a  very  thick  film 
was  observed  at  the  basal  face  (0001)  under  the 
same  conditions  (Fig.  6).  The  presence  of  air  signif¬ 
icantly  increased  the  liquid  film  thickness  at  the 
(0001)  face. 

NUCLEAR  MAGNETIC  RESONANCE 

Clifford  (1967),  using  Nuclear  Magnetic  Reso¬ 
nance  (NMR)  measurements  on  ice,  observed  a 
very  narrow  line  (-0.2  G)  that  he  attributed  to 
some  form  of  liquid  water  remaining  in  ice  down 
to  temperature  of  at  least  -12°C.  Kvlividze  et  al. 
(1970)  demonstrated  that  this  narrow  resonant 
line  is  distinct  from  the  NMR  lines  corresponding 
to  protons  in  water  and  bulk  ice  and  is  caused  by 
the  liquid-like  layer  on  the  ice  surface.  Afterwards, 
NMR  was  used  for  studying  the  properties  of  the 
liquid-like  layer  on  the  ice  surface  by  Bell  et  al. 
(1971),  Kvlividze  et  al.  (1974),  Barer  et  al.  (1977) 
and  Ocampo  and  Klinger  (1983a,b).  The  most 
comprehensive  information  about  the  physical 
properties  of  a  liquid-like  layer  was  obtained  by 
Mizuno  and  Hanafuza  (1987),  to  whom  we  are  go¬ 
ing  to  refer  most  often  for  this  reason. 

The  phenomenon  of  NMR  is  an  absorption  of 
electromagnetic  field  quanta  by  nuclei  that  pos¬ 
sess  a  magnetic  moment  and  are  placed  in  an  ex¬ 
ternal  magnetic  field  B  at  resonance  frequency  os0 

0)o  =  gB  (3) 

where  g  is  the  giromagnetic  constant. 

The  position  of  the  resonant  line,  its  shape  and 
half-width,  as  well  as  the  spin-lattice  relaxation 
time,  appear  to  be  quite  sensitive  to  the  properties 
of  a  medium  that  surrounds  the  nuclear  spins. 
Namely,  those  properties  are  the  average  distance 
and  relative  position  of  the  spins  in  space,  angular 
velocity  and  diffusion  coefficient  of  molecules  that 
contain  the  resonant  nuclei. 

Readers  who  are  interested  in  this  method 
should  see  the  marvelous  book  by  Abragam  (1961 ) 
and  also  the  paper  by  Burnett  and  Harmon  (1972), 

The  NMR  technique  applied  to  ice,  water  or  the 
quasi-liquid  layer  allows  us  first  to  distinguish 
three  resonant  lines  associated  with  absorption  by 
hydrogen  nuclei  in  those  three  media.  From  the 
intensity  of  the  liquid-like  layer  line,  we  can  judge 
the  approximate  thickness  of  a  particular  layer. 
Besides,  it  is  possible  to  determine  the  correlation 


time  of  rotational  motion  of  water  molecules  and 
a  self-diffusion  coefficient  (Mizuno  and  Hanafu- 
za  1987). 

Since  under  normal  conditions,  the  relative 
portion  of  molecules  in  the  liquid-like  layer  is 
negligibly  small,  it  is  necessary  first  to  increase 
the  "surface-to-volume"  ratio  of  the  ice  in  ques¬ 
tion  to  observe  NMR  from  this  layer.  The  easiest 
way  to  achieve  this  is  to  pulverize  water  at  a  low 
temperature.  As  a  result,  microscopic  ice  balls  are 
formed  with  a  diameter  about  100  pin  and  a  net 
surface  of  about  3-5x1 02  cm/g.  This  method  was 
used  for  instance  by  Mizuno  and  Hanafuza 
(1987).  Kvlividze  et  al.  (1974)  saturated  water 
with  air  prior  to  freezing.  This  leads  to  an  ap¬ 
pearance  of  a  great  number  of  microscopic  air 
bubbles  inside  the  ice,  with  a  net  surface  of  the 
order  of  ~104  cm2/g.  Frost  forming  on  walls  of  a 
Dewar  vessel  filled  with  liquid  nitrogen  has  an 
even  larger  specific  area  of  l-2xl05  cm2/g  (Kvliv¬ 
idze  et  al.  1970).  Placing  spherical  Teflon  particles 
(r  =  350  A)  (Kvlividze  et  al.  1974),  or  Aerosil 
particles  ( r  ~  100  A)  (Barer  et  al.  1977)  into  water 
before  freezing  produced  something  comparable 
to  the  "frost"  net  surface  of  ice  for  NMR  analysis. 
We  have  to  take  into  account,  of  course,  that  the 
structure  and  properties  of  the  liquid-like  layer 
on  tire  ice/ (solid  dielectric)  interface  might  be  dif¬ 
ferent  from  ice/ (water  vapor  +  air)  interface 
properties. 

Figure  7  shows  NMR  spectra  of  the  liquid-like 
layer  of  ice  measured  at  various  temperatures.  To 


Figure  7.  NMR  spectra  of  the  quasi-liquid  layer  of  ice 
and  water  observed  at  99.5  MHz.  The  signal  of  liquid 
water  was  taken  at  5°C.  Notice  the  difference  in  the  line 
width  between  liquid  water  and  the  quasi-liquid  layer 
(after  Mizuno  and  Hamfusa  1987). 


Figure  8.  Intensity  and  line  width  of  NMR  of  quasi¬ 
liquid  layer  against  inverse  temperature  (after  Mizuno 
and  Hanafusa  1987). 


increase  the  signal-to-noise  ratio,  the  spectra 
were  gathered  200  times  at  each  temperature.  The 
broad  signal  from  crystalline  ice  is  not  seen  with¬ 
in  the  range  of  the  20-kHz  observational  frequen¬ 
cy  (Mizuno  and  Hanafuza  1987). 

The  narrow  signal  was  never  detected  at  any 
temperature  if  bulk  ice  was  used  as  a  sample. 

As  is  clear  from  Figure  7,  the  line  width  and 
the  intensity  vary  with  temperature  and  their 
temperature  dependencies  are  shown  in  Figure  8. 
It  should  be  noted  that  the  line  width  of  die  liq¬ 
uid-like  layer  spectrum  at  -10°C  is  about  7  times 
that  of  ordinary  water  at  5°C. 

Since  the  temperature  dependencies  shown  in 
Figure  8  were  measured  as  temperature  was  in¬ 
creased  from  -100  to  -5°C,  and  since  at  T  ^  -10°C 
a  considerable  sintering  diminishes  the  total  ice 
surface,  the  actual  increase  of  the  NMR  signal  at  T 
>  -10°C  could  be  greater  than  the  one  shown  in 
Figure  8. 

The  temperature  dependence  of  the  correla¬ 
tion  time  for  rotational  motion  x  is  shown  in  Fig¬ 
ure  9.  This  time,  which  could  be  interpreted  as  the 
time  it  takes  a  water  molecule  to  be  reoriented  in 
a  liquid-like  layer  at  T  =  -10°C,  appears  to  be 
closer  to  that  time  for  water  (=10~12  seconds)  than 
for  ice  (2  x  10~5  seconds  [Jaccard  1959])  at  the  same 
temperature.  The  activation  energy  of  molecule 
reorientation  in  a  liquid-like  layer,  which  can  be 
determined  from  the  slope  of  the  curve  in  Figure 
9,  is  0.291  eV  and  is  also  quite  different  from  the 
corresponding  value  for  the  bulk  ice  (0.575  eV 
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Figure  9.  Correlation  time  for  rotational  motion 
X  of  quasi-liquid  layer  against  inverse  tempera¬ 
ture  ( after  Mizuno  and  Hanafusa  1987). 


[Jaccard  1959]).  Therefore,  comparing  the  rota¬ 
tions  of  a  water  molecule  in  water,  the  liquid-like 
layer  and  ice,  we  could  conclude  that,  according 
to  this  property,  the  liquid-like  layer  is  intermedi¬ 
ate  between  ice  and  water,  but  it  nevertheless 
resembles  water  more  closely. 

However,  if  we  try  to  compare  self-diffusion 
coefficients  in  water,  the  liquid-like  layer  and  ice 


(see  Fig.  10),  we  will  come  to  an  opposite  conclu¬ 
sion,  because  this  coefficient  at  0°C  in  the  liquid¬ 
like  layer  turns  out  to  be  =  3xl0~12  m2/s,  which  is 
closer  to  ice  («10-14  m2/s)  than  to  water  (2xl0-9 
m2/s).  At  the  same  time  the  activation  energy  of 
self-diffusion  in  the  liquid-like  layer  (0.24  eV)  is 
closer  to  that  of  water  (-0.24  eV)  than  to  ice  (0.65 
eV). 

Barer  et  al.  (1977)  studied  NMR  from  a  liquid¬ 
like  layer  forming  at  the  interface  between  ice  and 
Aerosil  microspheres,  i.e.,  on  a  frozen  aqueous 
suspension  of  Aerosil.  Like  Mizuno  and  Hanafu- 
za,  they  tracked  the  existence  of  the  liquid-like 
layer  much  further  down  on  the  temperature  scale 
than  measurements  that  could  be  obtained  from  X- 
rays  or  ellipsometry  (see  Fig.  11).  The  results  of 
their  estimations  of  the  liquid-like  layer  thickness 
at  various  temperatures  are  given  in  Table  1. 

Table  1.  Temperature  dependence  of  the 
thickness  of  a  liquid-like  layer  on  the  ice/ 

Aerosil  interface  (after  Barer  et  al.  1977). 

T{°  C)  -12  -10-8-6-4  -2 

1(A)  11  12  14  20  27  52 

hi  the  same  work  (Barer  et  al.  1977),  the  viscos¬ 
ity  of  a  liquid-like  layer  on  the  ice/quartz  interface 
was  measured.  At  T  =  -10°C  the  viscosity  ap¬ 
peared  to  be  q  =  0. 17  P,  which  exceeds  the  viscosity 
of  supercooled  water  at  the  same  temperature  by 
an  order  of  magnitude. 


Figure  10.  Diffusion  coefficient  of  quasi-liquid  layer  against  Figure  11.  Temperature  dependence  of  the  relative 
inverse  temperature  (after  Mizuno  and  Hanafusa  1 987).  concentration  of  the  mobile  aqueous  suspension  of 

Aerosil  (after  Barer  et  al.  1977). 
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ELECTRICAL  CONDUCTIVITY 
OF  THE  ICE  SURFACE 


As  we  have  seen  in  the  previous  sections,  the 
surface  layer  of  ice  differs  from  the  bulk  in  its 
structure,  optical  properties,  molecule  reorienta¬ 
tion  time  and  self-diffusion  coefficient.  It  is  there¬ 
fore  also  natural  to  expect  the  ice  surface  to  be  dif¬ 
ferent  from  the  bulk  in  its  electrical  properties. 

We  can  name  several  possible  things  that  could 
change  electrical  properties  of  ice  near  the  surface. 
The  first  is  an  increase  (or  decrease)  of  the  charge 
carrier  concentrations  within  the  screening  length, 
assuming  that  there  is  a  potential  difference  <p0 
between  tire  surface  and  the  bulk.  Such  a  differ¬ 
ence  may  arise  because  of  differences  in  the  chem¬ 
ical  potentials  of  charge  carriers  in  the  bulk  and  on 
the  surface  (electrochemical  potential  must  be  uni¬ 
form  everywhere)  or  because  of  adsorption  of  elec¬ 
trically  charged  impurities  that  compensate  the 
space  charge  of  ice  at  the  surface.  The  latter  arises 
merely  from  differences  in  solubilities  of  positive 
and  negative  impurity  ions  in  ice.  Second,  charge 
carrier  mobilities  and  the  mechanisms  of  interac¬ 
tion  between  them  might  be  different  in  the  disor¬ 
dered  surface  or  quasi-liquid  layer. 

In  fact,  the  ice  surface  has  an  extremely  high  net 
electrical  conductivity  that  was  discovered  back  in 
the  1950s  and  has  been  studied  intensely  over  the 
past  30  years.  At  temperatures  roughly  above 
-15°C  in  pure  ice,  the  d.c.  conductivity  of  the  sur¬ 
face  exceeds  that  of  the  bulk.  A  measurement 


ing  on  the  surface  from  the  bulk  currents.  Usual¬ 
ly,  this  is  done  by  means  of  so-called  guard  rings 
(Fig.  12).  A  guard  ring,  as  can  be  r  'en  from  this 
figure,  is  separated  from  the  central  electrode  by 
a  thin  insulating  ring.  The  guard  ring  and  the 
central  electrode  have  nearly  the  same  potential, 
which  provides _an  orientation  of  electrical  field 
strength  vector  E  (and  consequently  of  currents) 
parallel  to  the  surface  of  the  cylinder.  In  this  set¬ 
up  all  subsurface  currents  flow  through  the 
guard  ring,  and  their  magnitude  is  proportional 
to  the  voltage  drop  across  the  resistor  Rs.  In  this 
schematic,  Rs  and  R%  must  be  much  smaller  than 
tire  impedance  of  the  specimen  Zs  in  the  whole 
range  of  the  investigated  frequencies. 

A  great  majority  of  researchers  have  used  the 
configuration  shown  in  Figure  12,  although 
sometimes  square  ice  samples  were  used  instead 
of  cylindrical  ones.  Caranti  and  JQlingworth 
(1983a)  used  the  same  idea  (making  E  parallel  to 
the  surface  and  the  use  of  a  guard  ring),  although 
they  employed  a  semi-spherical  electrode  frozen 
into  ice. 

Jaccard  (1966, 1967)  suggested  that  four-point 
electrodes,  arranged  in  pahs  at  the  opposite  sides 
of  a  thin  ice  specimen,  be  used  to  determine  the 
surface  and  bulk  conductivities.  However,  we 
shall  refrain  from  discussing  this  here  because  of 
serious  criticism  of  this  method  expressed  by 
Turner  and  Stow  (1986), 


of  the  surface  conductivity  is  one  of  tire  sim¬ 
plest  and  most  accessible  ways  to  study  tire 
physical  properties  of  tire  ice  surface.  In  fact, 
a  simultaneous  measurement  of  tire  absolute 
value  of  the  surface  conductivity  a5  (with  a 
dimension  of  fl-1)  and  of  the  specific  volume 
electrical  conductivity  of  the  surface  layer  ass 
(with  a  dimension  of  Gr'm-1)  would  allow 
determination  of  tire  layer's  ‘hickness  L 

L  =  ■  (4) 

Measurement  of  charge  carrier  mobilities' 
values  p,  and  their  concentrations  nt  is  essen¬ 
tial  to  understanding  the  nature  and  the 
structure  of  tire  quasi-liquid  layer,  as  well  as 
for  determining  the  activation  energies  of  g, 
and  Researchers  hr  this  field  have  investi¬ 
gated  the  electrical  conductivity  and  dielec¬ 
tric  permittivity  of  the  surface  layer,  as  well 
as  the  Hall  effect. 

Tire  major  methodical  problem  in  these  ex¬ 
periments  is  the  separation  of  currents  flow- 


lo  power  supply 


insulating 

ring 


to  measuring  device 


guarding  ring 


Figure  1 2.  Typical  configuration  for 
measurements  of  surface  and  bulk 
conductivity. 
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Table  2.  Surfa.e  conductivity  of  ice  . 


T 

ra 

f 

(Hz) 

0, 

(O-i; 

Es 

(eV) 

Type  of  ice 

Reference 

-10 

d.c. 

3xl0-8 

0.84  ±0.1 

pure  polyc  ryslals 

Camp  et  al.  (1969) 

-10 

d.c. 

3.9x10  >° 

1.29 

pure  single  crystals 

Bullemer  and  Riehl  (1966) 

-10 

d.c. 

3.9x10-*° 

2  1.3 

3x10“*  MHF 

-10 

d.c, 

(1.5-6)xl0-*‘ 

1.52 

pure  single  crystals 

Maeno  (1973) 

-10 

d.c. 

(l-6)xl0-» 

1.19 

pure  single  crystals 

Maeno  and  Nishimura  (1978) 

-10 

d.c. 

6xl0-« 

0.44 

1.7x10“*  M  HH 

-10 

0.1 

4xl0-9 

1.0 

pure  polycrystals 

Caranti  and  Illingworth  (1983a) 

-10 

0.1 

6X10-6 

0.55 

10J  M  NaCl  polycrystals 

-10 

107 

2xl0-8 

0.78 

pure  poiycrystals 

Caranti  and  Lamfri  (1987) 

-10 

d.c. 

-  lO-5 

1.43 

pure  single  crystals 

Maidique  et  al.  (1971) 

-10 

d.c. 

1.5x10-*° 

1.29 

pure  single  crystals 

Ruepp  and  Kass  (1969) 

-11 

d.c. 

(1.4-6.3)xl0>° 

— 

pure  single  crystals 

Jaccard  (1967) 

-53 

d.c. 

-  1.5X10-*3' 

- 

polycrystalline  ice  film 

Chrzanowski  (1988) 

-10 

d.c. 

»  4x10  -8* 

deposited  onto  a  dielectric 
substrate 

*  Calculated  from  Chrzanowski's  data  assuming  that  the  electrode's  length  was  5  mm. 
t  Calculated  from  Chrzanowski's  data  assuming  that  the  activation  energy  is  1.2  eV. 


Even  a  small  amount  of  impurities  in  ice 
strongly  affects  the  surface  conductivity  as.  This 
accounts  for  the  considerable  scatter  in  the 
experimental  data  obtained  on  as.  Data  or.  low- 
frequency  and  d.c.  surface  conductivities  of  pure 
and  doped  ice  at  -10°C  are  summarized  in  Table 
2.  As  we  can  see  from  this  table,  os,  even  for  pure 
monccrvstallinp  ice,  assumes  values  in  a  wide 
range,  from  lG'a  Q~]  (Maeno  and  Nishimura 
1978)  up  to  ltr9!}-1  (Maidique  et  al.  1971).  Data 
on  activation  energies  Es  for  the  same  crystals 
vary  from  a  minimum  value  of  1.19  to  1.52  eV, 
with  an  average  vaiue  of  about  1.3  eV.  Even  a 
slight  doping  of  ice  leads  to  a  considerable  in¬ 
crease  in  crs  and,  as  a  rule,  a  decrease  of  £s. 


Typical  frequency  dependencies  of  ice  bulk 
and  surface  conductivities  are  plotted  in  Figure 
13.  This  plot  shows  that  the  frequency  depen¬ 
dence  of  the  bulk  conductivity  ag  follows  De¬ 
bye's  dispersion,  which  is  well  known  for  ice 

(Gm- an)  (xdU>)2  , 

an  ~  Co  +  -v - (5) 

1  +  (tdm)2 

where  o0  is  the  low-frequency  limit  of  the  con¬ 
ductivity  (or  d.c.  conductivity)  determined  by 
minority  charge  carriers.  In  pure  ice  at  -10°C, 
those  are  H3Q+  and  OH“  ions.  c„  is  the  high- 
frequency  limit  of  the  conductivity,  equal  to  the 
sum  of  all  partial  conductivities  of  the  charge 


O 

£ 

o  -  bulk  (crB) 

x  -  surface  (a,) 

•  -  surface  with  thick  frost 

■  -  surface  after  annealing. 


Figure  13.  Variations  of  the  real  part  of  conductivity  of 
pure  polycrystalline  ice  at  -10‘C  (after  Caranti  and 
Illingworth  1983). 
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Ogs  =  10-1  ft'1™-1. 


(7) 


carriers.  Under  the  conditions  considered,  <s„ 
practically  coincides  with  the  conductivity  of  the 
majority  charge  carriers — Bjerrum  defects,  that  is 
L-  and  D-defects.  co  is  the  angular  frequency  and 
ip  is  the  Debye  relaxation  time.  I  refer  readers 
who  are  more  interested  in  the  electrical  proper¬ 
ties  of  bulk  ice  to  the  latest  review  on  this  subject 
(Petrenko  1993).  Let  us  just  remember  that  a  small 
(compared  to  a„)  value  of  CTo  results  from  lattice 
polarization  arising  from  protonic  charge  carri¬ 
ers'  motion.  To  prevent  the  polarization  from  in¬ 
creasing  to  infinity,  ions  and  Bjerrum  defects 
have  to  move  in  a  consistent  manner  to  keep  their 
fluxes  equal.  An  absence  (or  reduction)  of  any 
noticeable  difference  between  o„  and  Cq  would 
imply  either  that  such  a  correlation  in  the  motion 
of  ions  and  Bjerrum  defects  is  absent  (or  reduced) 
or  that  partial  conductivity  of  ions  is  close  to  the 
conductivity  of  Bjerrum  defects. 

As  cait  be  seen  from  Figure  13,  the  frequency 
dependence  of  as  is  much  less  pronounced  than 
that  of  Ob.  This  is  especially  true  for  the  ice  sur¬ 
face  after  annealing,  i.e.,  the  smoothest  surface. 
Since  some  fraction  of  the  current  flowing  along 
the  surface  inevitably  penetrates  the  bulk,  the 
remaining  slight  frequency  dependence  can  be 
ascribed  to  these  accompanying  bulk  currents. 
The  rougher  the  surface  is,  the  greater  is  the 
"bulk"  contribution  of  the  "surface"  conductivi¬ 
ty,  in  excellent  conformity  with  the  results  shown 
in  Figure  13.  Thus,  we  can  tell  that,  at  least  in  the 
frequency  range  of  10-1  to  3xl04  Hz,  there  is  ei¬ 
ther  no  correlation  in  tire  motion  of  ions  and  Bjer¬ 
rum  defects  in  the  subsurface  layer  at  all  (as  is  the 
case  in  water)  or  that  this  correlation  is  weak. 

An  alternative  explanation  for  the  weak  fre¬ 
quency  dependence  of  cts  can  be  that  the  cor¬ 
responding  dielectric  relaxation  time  for  tire  sub¬ 
surface  layer  Tqs  might  be  several  orders  of  mag¬ 
nitude  more  than  that  of  the  bulk  To-  Therefore 
we  would  search  for  d  e  Debye  dispet  -ion  '->r  die 
liquid-like  layer  in  a  very-high-frequency  i„rve, 
v.-here  tire  standard  methods  for  measuring  a  are 
not  applicable. 

Let  us  try  to  estimate  the  specific  conductivity 
of  the  surface  layer  ass  (dimensions  of  ft-'m-1), 
proceeding  from  the  data  in  Table  2  and  the  re¬ 
sults  for  the  liquid-like  layer  thickness  L  given  in 
the  previous  sections. 

<*.  =  Sl.  .  (6) 

Thus,  taking  L  =  1.2xl0~9  m  (7’  =  -10C°  [Barer  et 
al.  1977])  and  as  =  10~10  ft-1,  we  obtain 


If  we  substitute  extreme  values  of  as  for  pure 
ice  from  10“n  ft-1  up  to  3x10^  ft-1,  we  shall  get 

Css l=  (lO^-SxlO1)  ft'1  nr1.  (8) 

At  the  same  temperature,  the  bulk  low-fre¬ 
quency  conductance  c0  falls  in  the  limits  (2.5x 
10-8  to  6.4xl(H0)  ft^irr1  (Petrenko  1993);  i.e.,  it  is 
from  6  to  11  orders  of  magnitude  less  than 
from  eq  8.  The  bulk  high-frequency  conductivity, 
c«  =  1 .6xl0~5  ft-1  m-1,  is  also  much  less  than  c^.  If 
we  assume  that  there  is  no  correlation  in  the  mo¬ 
tion  of  ions  and  Bjerrum  defects,  then  the  charac¬ 
teristic  dielectric  relaxation  time  of  the  surface 
layer  would  be  the  Maxwell  relaxation  time  Tms 

Tm s  =  -  (10-u  -  3x10-*)  s  (9) 

Os s 

which  corresponds  to  the  frequency  range 

/=— 1 — - (5xl07-  1.6xlOn)Hz  (10) 

271  IMS 

where  £o  is  the  vacuum  dielectric  pennittivity 
and  3.2. 

Unfortunately,  I  am  not  aware  of  any  attempts 
to  measure  os  as  a  function  of  frequency  in  the 
specified  range,  which  could  give  more  accurate 
values  of  in  conformity  with  eq  9.  If  in  the 
future  such  measurements  are  carried  out,  then 
this  will  also  allow  more  accurate  determination 
of  the  liquid-like  layer  thickness  L  using  eq  6. 

The  result  of  estimating  the  charge  carrier  con¬ 
centration  in  lire  liquid-like  layer  using  the  fol¬ 
lowing  formula  is  quite  interesting 

it,  =  SasL- .  (11) 

yaet 

If  we  take  (%  »  3xl0J  ft-'nr1  and  p  »  lO-7  m2  V-1 
s_1  for  ions  and  -  iu  m2  V-1  s_1  for  Bjerrum  de¬ 
fects  (Petrenko  1993),  then  ar  estimation  using  eq 
11  for  the  n  of  ions  and  Bjerrum  defects  gives 
3xl027  nr3  and  5xl028  nr3  respectively.  (We 
used  Cj  =  0.62e  and  e3  =  0.3Se  [Hubmann  1979].) 
These  are  extremely  large  values  if  we  recall  that 
the  concentration  of  water  molecules  in  ice  is 
about  3xl02S  m"3. 

Tie  above  estimations  can  be  much  more  accu¬ 
rate  if  we  employ  for  our  calculations  the  correla¬ 
tion  time  for  rotational  motion  of  water  mole¬ 
cules  at  the  ice  surface,  found  in  NMR  experi¬ 
ments  by  Mizuno  and  Hanafuza  (1987)  (see  Fig. 
9),  In  bulk  ice,  this  time  coincides  with  the  dielec- 


trie  relaxation  time  To  and  is  related  to  the  con¬ 
ductivity  0„  by  the  following  formula 

to  =  <I>  (- ica-  +  2a\=0  —  (12) 

lei  el)  el  el 

where  <P  =  3.85  k^T r(HV  r^,  =  2.76  A  is  oxygen- 
oxygen  distance,  "ion"  and  "B"  subscripts  corre¬ 
spond  to  ionic  and  Bjerrum  defect  conductivity 
and  t?!  and  e3  to  the  electrical  charges  of  ions  and 
Bjerrum  defects.  Substituting  To  a  10"'°  s  (from 
Pig.  8),  T  =  263  K  and  e3  =  0.38e  into  eq  12,  we 
obtain  for  the  net  conductivity  of  the  surface  layer 

or**  llil-' nr1  (13) 

i.e.,  a  value  close  to  the  upper  limit  found  in  eq  8. 
The  corresponding  concentration  of  Bjerrum  de¬ 
fects  in  the  surface  layer  according  to  eq  11  will  be 

«i3  =  2xl028  nr3.  (14) 

If  there  is  no  correlation  between  fluxes  of  ions 
and  Bjerrum  defects  in  tire  liquid-like  layer,  we 
should  use  eq  9.  Then 

Oss  =  0.3  SI"1  nr',  -  5xl026  in-3.  (15) 


Temperature  (UC) 

o  -5  -10  -15  -20  -25 


Figure  14.  Temperature  dependence  of  surface 
conductivity  of  the  (0001)  plane  for  pure  and 
doped  ice.  "O-lke"  stands  for  ice  from  O-ike  pond 
near  the  Syowa station  in  Antarctica,  "Mendenhall" 
stands  for  ice  front  Mendenhall  Glacier,  Alaska  (after 
Maenoand  Nishimura  1978). 


Thus,  we  can  conclude  that  the  charge  carrier 
concentration  in  the  ice  surface  exceeds  by  many 
orders  of  magnitude  that  of  the  bulk  and  in  some 
cases  might  correspond  to  a  superionic  state,  i.e., 
a  state  in  which  all  (or  nearly  all)  molecules  are 
ionized  and  mobile. 

Caranti  and  Illingworth  (1983a)  remark  in 
their  paper  that,  in  a.c.  measurements  of  real  and 
imaginary  parts  of  the  complex  surface  conduc¬ 
tivity,  the  phase  angle  was  always  small,  show¬ 
ing  surface  currents  to  be  predominantly  resis¬ 
tive.  This  is  exactly  what  must  be  the  case  for  an 
ordinary  conductivity  (without  correlation  be¬ 
tween  charge  carrier  fluxes)  at  frequencies  much 
less  than  td"'. 

Temperature  dependencies  of  the  surface  con¬ 
ductivities  of  pure  and  doped  ice  are  plotted  in 
Figures  14  and  15.  The  behavior  of  this  depen¬ 
dence  is  different  at  temperatures  above  and  be¬ 
low  -6°C.  It  is  quite  possible  that  a  sharp  rise  of  as 
at  T  £  -6°C  is  associated  with  an  appearance  of 
the  liquid-like  layer  on  the  ice  surface,  which  is 
detected  by  optical  ellipsometry  and  X-ray  dif¬ 
fraction  methods  (see  previous  sections).  An 
interesting  observation  is  the  difference  in  as  cor¬ 
responding  to  different  crystallographic  surfaces 
of  ice  (Fig.  15).  The  conductivity  of  the  basal 

Temperature  (°C) 

0  -5  -10  -15  -20 


Reciprocal  Temperature  (10  -'K  _l) 

Figure  15.  Surface  conductivity  of  ice  as  a  func¬ 
tion  of  crystallographic  orientation  and  tempera¬ 
ture.  Ice  single  crystals  were  taken  from  the 
Mendenhall  Glacier,  Alaska  (after  Maeno  1973). 
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Figurel6.  Variations  of  surf  ace  conductance  after  evacu¬ 
ation  of  pure  water  for  "pure"  natural  ice  (O-ikcpond). 
Dotted  lines  indicate  expected  variations  of  the  conduc¬ 
tance  calculated  from  measured  temperature  of  the  ice 
surface  (after  Maeno  and  Nishimura  1978). 

(0001)  surface  was  maximum.  This  result  is  in 
agreement  with  optical  measurements  by  Elbaum 
et  al.  (1993)  but  is  opposed  to  the  results  of  optical 
ellipsometry  (where  liquid-Uke  layer  thickness  is 
greater  on  the  prismatic  (1010)  surface)  and  X-ray 
diffraction  (where  the  difference  between  the 
(0001)  and  (1010)  surface-,  is  negligible). 

Maeno  (1973)  and  Maeno  and  Nishimura  (1978) 
found  that  an  intense  pumping  out  of  water 
vapors  from  above  the  ice  surface  changes  the 
surface  conductivity  considerably  (Fig.  16).  Pump¬ 
ing  out  of  water  vapors  accelerates  evaporation 
from  the  surface  and  its  cooling.  Aside  from  this 
trivial  effect,  a  lower:  j  V  water  vapor  prcsaui : 
above  the  ice  increases  ti.e  surface  conductivity. 
The  phenomenon  discovered  by  Maeno  shows  the 
degree  to  which  the  conductivity  of  the  surface 
layer  is  sensitive  to  physical  conditions  at  the  ice 
interface. 

An  interesting  result  was  obtained  in  the  work 
by  Chrzanowski  (1988),  who  investigated  the  d.c. 
conductivity  of  thin  ice  films  deposited  onto  a 
high  resistance  dielectric  substrate  (high  quality 
glass  or  quartz  plate).  Figure  17  demonstrates  that 
the  conductance  of  an  ice  film  grows  rapidly  with 
the  film's  thickness  until  it  thickens  to  0.3-0.5  pm, 
Then,  the  rate  of  tire  conductance  growth  dimin¬ 
ishes  considerably.  This  result  implies  that  the  lay- 
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Current 


pm 

Figure  17.  Typical  thickness  dependence  of  the  mea¬ 
sured  current  intensity  obtained  for  ice  films  deposited 
onto  a  glass  substrate  of  temperature  220  K.  Three 
different  experimental  runsureshown  (after  Chrzanowski 
1988). 

er  close  to  the  ice /dielectric  interface  has  a  much 
larger  specific  conductivity  than  the  bulk  and  that 
the  thickness  of  this  layer  is  about  0.3-0.5  pm.  This 
thickness  (taking  also  into  account  the  temper¬ 
ature  of  220  K)  exceeds  considerably  values  found 
in  other  works.  I  be<ieve  that  this  quantity  can  be 
associated  only  with  the  screening  length  in  ice 
which  was  considered  in  detail  by  Petrenko  and 
Ryzhkin  (1984). 

Caranti  and  Lamfri  (1987)  made  a  very  interest¬ 
ing  attempt  to  determine  charge  carrier  mobilities 
in  the  ice  surface  layer.  The  authors  used  for  this 
purpose  measurements  of  the  Hall  effect  and  sur¬ 
face  conductivity.  Ovring  to  the  very  small  mo¬ 
bilities  of  charge  carriers  in  ice,  the ;  asuremerrts 
of  the  Hall  euect  are  quite  complicated.  Caranti 
and  Lamfri  used  an  a.c.  setup  with  two  distinct 
frequencies  of  the  magnetic  field  (50  Hz)  and  the 
driving  current  (107  Hz),  measuring  the  Hall  effect 
at  tire  total  frequency  157  Hz.  Figure  18  shows  the 
temperature  dependence  of  tire  carrier  mobilities 
in  the  ice  surface  determined  by  Carai  *i  and  Lam¬ 
fri.  An  unusually  large  value  of  mobility  (3xl0~* 
m2  V-1  s_1)  is  striking.  This  is  three  orders  of  mag¬ 
nitude  greater  than  tire  mobility  of  ions  in  water 
and  ice  and  four  orders  greater  than  the  mobility 
of  Bjerrum  defects  in  ice  (Petrenko  1993).  More¬ 
over,  such  high  mobility  can  not  be  attributed  to 
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Figure  18.  Hall  mobility  g  (m2  V'1  s'1)  at  157  Hz  as  a 
function  of  temperature  for  normal  and  rubbed  sur¬ 
faces.  The  magnetic  field  By  was  0.38  T  (after  Carantiand 
Lamfri'1987). 


Nonetheless,  the  result  of  Caranti  and  Zamfri 
is  indirectly  confirmed  by  earlier  measurements 
of  Bullemer  and  Riehl  (1968),  who  tried  to  mea¬ 
sure  a  bulk  Hall  effect  in  ice  and  obtained  p  = 
1.4xKHm2V_1s_!  at  T  =  -8°C.  Although  Bullem¬ 
er  and  Riehl  used  guard  rings  at  the  Hall  elec¬ 
trodes  to  eliminate  the  effect  of  surface  currents, 
they  apparently  did  not  succeed  because  they  did 
observe  the  Hall  effect,  which  is  forbidden  for 
bulk  ice  in  principle  (Sokoloff  1973,  Gosar  1974). 
Hence,  Bullemer  and  Riehl  also  determ  ined  tire 
surface  Hall  mobility.  It  is  remarkable  that  their 
result  (although  they  used  d.c.  measurements)  is 
in  a  good  agreement  with  the  Caranti  and  Lamfri 
result.  Then,  with  what  is  this  extremely  high 
mobility  associated?  Perhaps  with  electrons  in 
the  ice  surface?  We  do  not  know  the  answer. 

Finally,  Camp  and  Halchin  (1979)  studied  the 
electrical  properties  of  ice/ice  (grain  boundaries) 
and  ice/quartz  interfaces  in  the  temperature 
range  from  0  to  -20°C  and  in  the  frequency  range 
from  d.c.  to  107  Hz.  They  found  that  the  widths  of 
the  ''electrical"  interfaces  (ice/ice  and  ice/dielec¬ 
tric),  estimated  from  the  conductivity  mea¬ 
surements,  depend  on  the  frequency.  The  higher 
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Table  3.  Summary  of  estimates  of  grain  boundary  thickness  8  (after  Camp 
and  Halchin  1979).  The  measuring  frequency  is  indicated  in  parentheses. 


Sample _ Type  of  measurement  Effective  tfrain  boundary  thickness 


Pure  no.  2 

From  high  frequency 

-rc 

6  (5  x  1U6)  100  A  (order  of  magnitude) 

Pure  no.  3 

From  capacitance  data 

8  (K)6)  =  700  A 

Pure  no.  1 

From  d.c.  compared  to  doped  ice 
From  capacitance  data 

8  (d.c.)  =  2400  A  (order  of  magnitude) 

8  (II)6)  -  ltXH)  A 

Doped 

From  capacitance  data 

-rc  -i4°c 

8(106)  =  12pm  3  pm 

8  (107)  =  5  pm  0.3  pm 

ions  at  all.  Even  in  tb~  '  i  of  an  activation-free 
diffusion  motion  c/  >  ms,  their  mob'lity  cannot 
exceed  the  value 

(l  •=  (16) 

k\}F 

where  v  -  frequency  of  ion's  hops  between  mol¬ 
ecules 
e  = charge 

a  =  length  of  a  hop  (=  10~10  m). 

To  explain  p  =  3  x  10-4  m2  V-1  s-1,  we  have  to 
accept  that  v  =  1015  s-1,  which  is  much  greater 
than  all  possible  phonon  frequencies  in  ice  (as 
roll  as  in  any  other  material). 


the  frequency,  the  thinner  the  surface  layer  in 
which  the  dielectric  relaxation  time  is  small 
enough  for  water  molecules  to  be  able  to  reorient. 
That  is,  the  surface  layer  is  not  homogeneous:  the 
charge  carrier  concentration  is  greater  near  tire 
surface  and  decreases  away  from  the  surface  into 
the  bulk  (see  Table  3,  which  is  a  summary  of 
results  of  Camp  and  Halchin). 

SURFACE  CHARGE  AND 
SURFACE  POTENTIAL 

The  study  of  the  surface  potential  of  ice  and 
electric  charges  at  tire  ice  surface  has  quite  a  long 
history  and  continues  to  attract  the  attention  of 
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numerous  scholars.  Such  interest  in  the  electrostat¬ 
ic  properties  of  the  surface  of  ice  is  stimulated  by 
both  fundamental  problems  associated  with  it  and 
certain  practical  issues.  Fundamental  questions, 
which  the  study  of  the  surface  potential  and  the  ice 
surface  might  help  to  resolve,  include  the  issue  of 
the  structure  of  the  surface  of  ice.  For  instance,  if 
all  the  molecules  at  the  ice  surface  are  oriented 
"proton-outwards,"  as  some  theoretical  models 
suggest  (Fletcher  1968),  this  must  lead  to  a  positive 
charge  at  the  surface  X*  and  a  positive  potential  of 
the  ice  surface  <ps.  Both  X*  and  <ps  are  determined  by 
the  thickness  of  the  surface  layer  d  and  the  fraction 
of  the  oriented  water  molecules  contained  in  this 
layer.  Thus,  measurement  of  X,  and  <ps  in  conjunc¬ 
tion  with  other  data  (d,  for  example)  might  pro¬ 
vide  valuable  information  on  the  microscopic 
structure  and  nature  of  the  ice  surface  layer. 

The  practical  issues  stimulating  the  study  of 
electrical  properties  of  tine  ice  surface,  incorporate, 
first  of  all,  the  problem  of  atmospheric  electricity 
(see  details  in  Hobbs  1974)  as  well  as  the  problem 
of  adhesion  and  friction  of  ice  and  snow  (Petrenko 
1994). 

Before  we  proceed  to  a  review  of  experimental 
data,  let  us  consider  basic  theoretical  relationships 
defining  the  electrostatic  properties  of  the  ice  sur¬ 
face  and  the  most  frequently  employed  experi¬ 
mental  techniques  used  to  determine  of  X^  and  <ps. 

Figure  19  shows  schematically  the  distribution 
of  the  electric  charge  density  p  in  ice  near  the 
surface  (x  =  0).  The  coordinate  scale  is  chosen  de¬ 
liberately.  First,  since  ice  is  a  conductor,  in  equilib¬ 
rium,  the  electric  field  within  the  bulk  must  be 
zero,  i.e.,  £(+  <*>)  =  0.  If  we  use  the  formula  for  the 


electric  field  strength  above  a  charged  surface 


£  *  is*. 

2e0 

we  shall  obtain  an  apparent  result 

r+“ 


£(+<*>) 


2eo 


p (x)dx  -  0. 


(18) 


(19) 


In  other  words,  the  sum  of  all  charges  in  the  sur¬ 
face  layer  of  ice  must  be  zero. 

A  very  thin  layer  of  the  surface  charge  per  se, 
with  the  density  Xso,  is  composed  of  ions  adsorbed 
by  the  surface  from  the  air,  electrons  entrapped  in 
the  surface  layer,  etc.  This  layer  is  screened  on  the 
inner  side  by  a  much  thicker  layer  of  a  space 
charge,  composed  of  mobile  charge  carriers,  i.e., 
predominantly  Bjerrum  defects  and  ions.  Using 
Poisson's  equation,  we  obtain 


dE  =  PW_ 
dx  £oe~ 


(20) 


and  as  dtp/dx  =  -£,  we  can  easily  find  the  relation¬ 
ships  between  p,  £  and  <p 

EW  =  I!  d%;  <p«  -  -  •  (2D 

JO  EoEoo  JO 

£«,  is  high-frequency  dielectric  permittivity  of 
ice, 

To  estimate  a  fraction  of  water  molecules  orient¬ 
ed  along  the  electric  field,  i.e.,  with  dipole  mo¬ 
ments  oriented  normally  to  the  surface,  we  have  to 
find  the  component  of  the  polarization  vector  p 


Figure  19.  Space  distribution  of  electric  charges  in  ice  near  the  surface. 
Scales  are  arbitrary. 


that  is  determined  by  the  reorientation  of  water 
molecules  (Pm) 


(22) 


where  «h20  a^d  n±  are  water  molecule  concentra¬ 
tion  in  ice  and  the  concentration  of  water  mole¬ 
cules  oriented  along  the  field  respectively.  Since 


/5M=£o(es-e«)£ 


■  (Hr) 


P&M  (23) 


then 


nj,  _  Pm  _  (  £s  ~ £- 
MHjO  Mnino  (e^MnujO 


P &)d%. 


(24) 


Here,_Ri  is  the  static  dielectric  permittivity^of 
ice  and  M  is  the  molecular  dipole  moment  (M  = 
8.66X10"30  °C  m  [Hobbs  1974]). 

This  molecular  "ordering"  reaches  its  maxi¬ 
mum  at  the  surface  where 


somehow  to  zero  when  x—> «  since  tire  ice  bulk  is 
disordered.  Then  the  surface  charge  Au,  associated 
with  this  polarization  is 

A.s  =  <W-dx  =  Px(  0).  (27) 

Jo  dx 

To  find  the  other  quantity  that  can  be  mea¬ 
sured  from  an  experiment,  namely,  the  surface 
potential  <(>s  (using  eq  21),  we  have  to  know  how 
tire  charge  p(x)  is  distributed  in  ice;  in  other 
words,  we  have  to  know  the  laws  of  screening  of 
an  electric  charge  in  ice.  The  screening  of  the 
charge  A,  will  occur  differently,  depending  on  the 
ratio  of  the  average  thermal  energy  kuT  to  the 
electrostatic  energy  of  the  charge  carriers  I  e,  <ps  I , 
If  the  surface  potential  is  small,  so  that 


kuT  »  I  ej  <ps  I  (28) 

then  in  this  case,  as  was  shown  by  Petrenko  and 
Ryzhkin  (1984) 


til  _  £s  “  1 

- * — - — — 

UH20  e«JVl  HH20 


(25)  £(x)  =  exp(-aeix)  +  exp(-aez;c)  (29) 

2eoe«  2eoes 


ksa  = 


P  (W- 


JO 


(26) 


Note  that  the  fraction  of  ordered  molecules  in 
eq  25  depends  only  upon  the  density  of  charge  L,a 
adsorbed  by  the  ice  surface.  In  other  words,  this 
result  will  still  hold  if,  owing  to  the  change  of 
temperature  or  doping,  the  concentration  of 
charge  carriers  changes,  resulting  in  a  change  of 
the  screening  charge  p(x)  in  the  bulk.  Therefore, 
eq  25  provides  a  remarkable  and  simple  way  to 
study  tire  ordering  of  water  molecules  in  the  sur¬ 
face  layer  through  a  measurement  of  tire  electric 
charge  A^,  adsorbed  in  a  thin  surface  layer.  The 
thickness  of  this  layer  is  much  less  than  the 
screening  radius  ar1. 

Above  we  have  considered  the  eiectric  charge 
absorbed  by  the  ice  surface  as  a  reason  for  the 
ordering  orientation  of  water  molecules.  There 
can  be  an  inverse  phenomenon  when  an  ordering 
of  water  molecules  on  the  ice  surface  generates  a 
surface  charge.  In  its  turn,  the  ordering  of  dipole 
moments  of  water  molecules  may  be  caused,  for 
example,  by  dipole-dipole  and  dipole-quadru- 
pole  interactions  (Fletcher  1968).  Suppose  this  or¬ 
dering  is  associated  with  polarization  P(x)  that 
reaches  its  maximum  at  tlie  surface  and  decays 


where  ee'J  and  as~2  are  tlie  two  screening  lengths 
determined  by  the  major  charge  carriers  (Bjerrum 
defects)  and  minor  charge  carriers  (ions)  respec¬ 
tively 


^-1  _  /  eoe-^uT 

V  Cs(»3+  >u) 

(30) 

a r2‘«J  - 

V  e^ini  +  m) 

(31) 

0> »  3.85  kvT  mo. 

(32) 

Since  «  Eg,  a  large  fraction  of  electric  field 
(and  potential)  decays  according  to  eq  29  at  the 
distance  ae  j  «  ae  2  •  This  happens  because  of 
screening  by  major  charge  carriers.  The  rest  of  the 
potential  decays  at  the  characteristic  length  a;  2 , 
determined  by  minor  charge  carriers. 

hi  the  case  of  strong  fields  (large  densities  of 
the  surface  charge),  when 

kuT  «  I  e,<ps  I  (33) 

the  screening  is  described  by  a  system  of  nonlin¬ 
ear  differential  equations,  which  apparently  does 
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not  have  an  analytical  sulution.  Petrenko  and 
Maeno  (1987)  have  found  a  solution  by  making 
an  approximation  that  ignores  the  difference  be¬ 
tween  fractional  values  of  the  electric  charges  of 
Bjerrum  defects  and  ions.  Assuming  £3  =  e\  =  q, 
we  obtained 

£  =  Jhm - 1 -  (34) 

2£ofc>«4  ^  -f  ^ 

where  the  screening  length  ae-1  now  depends  on 

^■sa 

-i  4ki}  Teo  £«*  /*3c\ 

m  = - .  (3b) 

q  Asa 

Equation  33  can  be  used  in  the  subsurface  region 
to  a  depth  where  tire  potential  I  <ps  £1 1  becomes 
comparable  to  k$T.  Further  down,  we  should  use 
eq  29. 

Before  proceeding  to  the  discussion  of  experi¬ 
mental  data,  I  shall  point  out  an  apparent 
relationship  that  exists  between  and  the  com¬ 
ponent  of  the  surface  conductivity  as,  which  is 
determined  by  the  layer  Xs  of  the  screening 
charge  of  the  protonic  carriers.  If  in  the  bulk  ( x  = 
+«>),  the  charge  carrier  concentrations  are  denot¬ 
ed  as  «j  (+«>),  then  the  corresponding  values  of 
concentrations  at  the  surface  will  be 

ni(0)a«,(Hexp(-g).  (36) 

In  reality  one  of  tire  four  types  of  drarge  carri¬ 
ers  always  dom mates,  hence  the  high-frequency 
surface  conductivity  is  given  by 

as  =  2  Id  Hi  I  tudx  =|A,Jpi .  (37) 

i=l  JO 

Petrenko  (1994)  and  Petrenko  and  Colbeck  (in 
press)  studied  an  electrification  arising  from  fric¬ 
tion  between  metal  and  dielectric  sliders  on  ice. 
Thin  films  rubbing  against  the  surface  of  a  rotat¬ 
ing  ice  cylinder  were  used  as  sliders.  The  experi¬ 
ments  were  conducted  in  the  -4.5  to-31.5°C  tem¬ 
perature  range  and  in  the  0.5-  to  8-m  s-1  velocity 
range;  tire  maximum  pressure  exerted  on  the  ice 
by  a  slider  was  1  kPa.  In  these  experiments,  we 
found  that  during  sliding,  tire  sliders  pick  up  a 
significant  positive  charge  from  ice.  In  tire  case  of 
metal  siiders  (stainless  steel  and  aluminum),  a 
large  potential  difference  of  1600  V  arises  be¬ 
tween  the  ice  and  tire  slider  (velocity  =  8  m/ s,  T  = 
-31.5°C).  When  dielectric  sliders  were  used  (un¬ 
der  the  same  conditions)  strong  electric  fields  £  = 


2.1xl0b  V  m_1  were  produced  at  the  ice/slider 
interface.  From  the  meas-urements  of  tire  current  1 
passing  through  an  electrometer  used  as  a  current 
meter,  we  found  that  a  coupled  ice  and  slider 
acted  as  a  charge  generator.  Thus 

/ =  —  =  Xs  —  (38) 

dt  ds 

where  Q  =  electric  charge 

X  -  surface  charge  density  which  the  slid¬ 
er  picks  up  from  the  ice  surface 
s  =  swept  area 
t  =  time. 

Similar  results  were  obtained  in  the  course  of  field 
experiments  on  sliding  of  alpine  skis  on  ice.  The 
value  of  tire  electric  charge  density  found  is  enor¬ 
mous:  Xs  =  1.6 xlO'2  C  nr'2.  But  even  this  value 
must  be  considered  as  the  low  limit  for  Xg  since, 
first,  the  slider  is  unlikely  to  pick  up  the  entire 
charge;  second,  if  the  slider  penetrated  to  sonre 
depth  into  the  ice,  some  fraction  of  the  unwanted 
screening  charge  can  be  entrapped  along  with  Xg. 

It  is  interesting  to  make  some  estimates  based 
on  the  value  of  Xs  that  we  found.  First,  such  a  high 
density  of  the  surface  charge  corresponds  to  0.01 
electron  charges  per  water  molecule  in  the  ice  sur¬ 
face.  Second,  let  us  estimate  the  fraction  of  water 
molecules  on  the  surface  that  are  directed  per¬ 
pendicularly  to  the  surface 

(39) 

"1120  M 

Substituting  X,  -  1.6xl0~z  C  m-2  into  eq  27,  we 
find  from  eq  39  that 


-4k- =0.2. 

«H20 


(40) 


And  remember  this  is  the  lower  limit!  It  means 
that  a  significant  fraction  of  water  molecules  at  ice 
surface  must  be  polarized. 

Finally,  using  eq  37  to  calculate  the  surface  con¬ 
ductivity  <js,  we  obtain 


as=  1.6xl0-2Cm~21.8xl0'8m2V-1s_1 


«  3xlO~10 £2_1  (41) 

where  the  value  of  the  mobility  of  L-defects 
(which  must  be  the  major  charge  carriers  screen¬ 
ing  a  positively  charged  surface)  at  T  =  -10°C  (Jac- 
card  1959)  was  used.  The  estimated  value  of  as 
found  in  eq  41  agrees  in  the  order  of  magnitude  to 
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experimental  data  at  T  =  -10°C  (see  Fig.  14  and  15). 
Although  it  seems  unlikely  that  in  such  a  way  all 
properties  of  cs  can  be  explained  (the  temperature 
dependence  of  (i4  and  os  is  different,  for  instance), 
the  conformity  found  is  hardly  purely  incidental 
and  implies  that  this  process  significantly  contrib¬ 
utes  to  the  surface  conductivity. 

Buser  and  Jaccard  (1978)  studied  a  charge  ex¬ 
change  between  very  small  (10-pm)  ice  particles 
and  various  metal  targets,  doing  the  experiments 
at  T  =  -45°C  and  collision  velocities  of  10  m/s. 
They  found  that  the  amount  and  the  sign  of  the 
charge  transferred  during  a  collision  depended 
upon  tire  work  function  of  the  metal  target  as  well 
as  the  value  and  the  sign  of  the  external  electric 
field,  in  which  tire  particles  were  polarized  prior  to 
collision.  To  interpret  their  experimental  results, 
Buser  and  Jaccard  assumed  the  presence  of  elec¬ 
tronic  states  at  the  ice  surface  in  an  analogy  to 
semiconductors.  These  states  participate  in  the 
electron  exchange  between  ice  and  metal  during  a 
short  (10-7  seconds,  in  the  described  experiments) 
collision.  The  lower  limit  of  the  density  of  electron 
surface  states  was  found  to  be  Xsa/e  2  1016  nr2.  If 
we  assume  that  the  amount  of  charge  AQ  that 
passes  between  the  ice  and  the  metal  during  the 
collision  is  proportional  to  the  difference  in  their 
work  functions  (Wj  -  Wrn),  then  TV;  can  be  defined 
as  the  point  of  intersection  with  the  horizontal  axis 
on  the  graph  of  AQ  versus  Wm.  In  this  manner, 
Buser  and  Jaccard  found  that  VVj  =  4.3  eV.  A  simi¬ 
lar  result  was  obtained  in  experiments  of  the  same 
kind  by  Buser  and  Aufdermaur  (1977).  Caranti  et 
al.  (1985),  who  used  the  same  method  to  de¬ 
termine  the  work  function  of  ice,  found  a  some¬ 
what  smaller  value  for  Wj  of  4.11±0.2  eV. 

Another  method  used  by  many  authors  is 
Kelvin's  method.  Figure  20  illustrates  the  princi¬ 
ple  of  measurement  of  the  difference  between  a 
work  function  of  the  studied  material  (in  our  case 
Wj,  the  work  function  of  ice)  and  that  of  an  oscillat¬ 
ing  electrode  Wm.  In  this  method,  ice  (or  another 
material)  is  connected  electrically  to  the  electrode 
that  is  oscillating  above  its  surface.  If  Wm  £  W,, 
electrons  flow  from  the  ice  into  the  metal  if  Wm  > 


W,  or  in  the  opposite  direction  if  Wm  <  Wi(  until 
the  equilibrium  potential  difference  Acp  is  estab¬ 
lished 

A(p  e  =  Wm  -  Wj .  (42) 

This  potential  difference  produces  an  electric 
field  in  the  space  between  the  ice  and  the  oscillat¬ 
ing  electrode,  which  forms  a  flat  capacitor.  When 
tire  electrode  oscillates,  the  capacitance  also  os¬ 
cillates,  resulting  in  an  alternating  electric  charge 
at  the  capacitor  and  an  electric  current  in  the  cir¬ 
cuit.  At  the  voltage  V  =  -A(p,  the  magnitude  of  the 
current  is  minimal,  which  allows  Wm  -  W;  to  be 
determined.  To  find  the  absolute  value  of  Wj,  we 
need  to  find  Wm  using  some  other  method  (from 
the  extrinsic  photo-efiect,  which  has  been  done 
already  for  all  pure  metals,  for  instance).  Kelvin's 
method  was  used  to  determine  the  work  function 
of  ice  by  Takahashi  (1969a,b,  1970),  Mazzeda  et  al. 
(1976),  Caranti  and  Illingworth  (1980,  1983a,b,c) 
and  Caranti  et  al.  (1985). 

The  approximate  value  of  Wj  found  in  these 
experiments  appeared  to  be  quite  sensitive  to 
such  ambient  conditions  as  humidity,  tempera¬ 
ture  and,  to  a  smaller  degree,  ice  doping.  A  char¬ 
acteristic  maximum  value  by  which  W;  was  al¬ 
tered  due  to  a  change  in  the  ambient  conditions 
was  1 8W;  I  <  200  meV.  Thus,  for  example,  Taka¬ 
hashi  (19/  w  observed  a  buildup  of  a  positive  po¬ 
tential  (30-80  mV)  at  the  ice  surface  at  the  begin¬ 
ning  of  the  condensation  process.  During  further 
condensation,  the  surface  potential  changed  its 
sign  and  reached  a  value  of  -250  mV.  The  effect 
was  more  pronounced  at  the  basal  (0001)  plane  as 
compared  to  the  prismatic  plane.  Condensation 
and  evaporation  at  the  ice  surface  were  produced 
through  its  heating  or  cooling,  respectively;  the 
ambient  temperature  was  maintained  constant. 

The  same  author  observed  a  negative  potential 
buildup  at  the  ice  surface  after  its  rubbing  (Taka¬ 
hashi  1969b).  This  result  is  in  qualitative  agree¬ 
ment  with  the  observations  of  Petrenko  and  Col- 
beck  (in  press),  although  we  found  a  much  larger 
effect. 

Mazzeda  et  al.  (1976)  studied  variance  in  Wj  as 
a  function  of  temperature.  In  these  experiments, 
tire  ice  surface  was  in  equilibrium  with  the  satu¬ 
rated  vapor  throughout  the  entire  temperature 
range  of  223-273  K;  thus,  the  effect  of  condensa¬ 
tion  and  evaporation  should  not  have  manifest¬ 
ed.  Figure  21  shows  a  typical  temperature  depen¬ 
dence  Of  AWjjy, 
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Figure  21.  Temperature  dependence  of  AW  im  =  l/|e| 
(W  m  -  W,)  (after  Mazzeda  etal.  1976). 
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Figure  22.  Change  in  ice  surface  potential  when  the 
surface  was  rimed  with  supercooled  droplets  (after 
Caranti  and  Illingworth  1983c). 
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The  arrow  in  F;p"re  21  indicates  the  tempera¬ 
ture  (=  -30°C)  at  v.n.Jt  there  is  a  transition  froma 
weak  linear  dependence  to  a  stronger  logarithmic 
dependence.  The  surface  potential  grew  more 
positive  relative  to  the  bulk  and  reference  elec¬ 
trode  when  it  approached  the  melting  point, 
which  can  be  explained  in  terms  of  the  formation 
of  a  double  charge  layer  near  the  surface,  as 
shown  in  Figure  19.  Such  charge  distribution  aris¬ 
es  if  the  molecules  at  the  surface  are  oriented 
oxygen-outwards,  entrapping  a  positive  com¬ 
pensating  charge  from  the  atmosphere. 

Caranti  and  Illingworth  (19  50, 19d3a,b,c)  stud¬ 
ied  variations  in  tire  surface  potential  when  the 


ice  surface  was  rimed  with  supercooled  droplets 
from  an  ultrasonic  generator  (of  50-gm  radius;  1- 
to  2-g  rh"3  water  content).  A  drastic  change  in  the 
surface  potential  was  observed  (Fig.  22).  The  ef¬ 
fect  of  doping  by  NH3,  NaCl  and  HF  (up  to  con¬ 
centrations  of  10~2  M)  on  the  amount  of  variation 
in  potential  was  insignificant.  The  authors  came 
to  the  conclusion  that  the  variation  in  the  surface 
potential  of  ice  during  its  riming  is  a  result  of  the 
disordered  ice  produced  by  rapid  freezing  of  the 
supercooled  water,  since  this  change  is  propor¬ 
tional  to  the  rate  of  freezing  and  relaxes  with  a 
time  constant  comparable  to  that  of  a  mechanical 
relaxation. 

SURFACE  ELECTROCHEMICAL 
POTENTIAL 

When  the  impurity  concentration  in  ice  varies 
from  place  to  place,  a  static  electric  field  arises 
from  inhomogeneity  in  the  distribution  of  impu¬ 
rities,  which  are  the  effective  sources  of  protonic 
charge  carriers  in  ice.  As  a  result  of  this  inhomo¬ 
geneity  of  impurity  concentration,  a  gradient  of 
the  charge  carriers'  concentrations  «,■  arises.  The 
fluxes  of  ions  and  Bjermm  defects  are  correlated 
since  they  alter  the  orientation  of  the  water  mole¬ 
cule  when  they  pass  through  it.  The  mathemati¬ 
cal  description  of  the  general  case,  when  both  dif¬ 
fusive  and  drift  fluxes  of  all  four  types  of  charge 
carriers  must  be  taken  into  account,  is  quite  com¬ 
plicated,  Nevertheless,  we  can  easily  make  rough 
estimations  of  the  magnitude  of  the  arising  elec¬ 
tric  fields,  using  a  simplified  consideration. 

Let  us  regard  the  case  of  on*1  *ype  of  charge 
carrier  of  electric  charge  e,-  and  concentration  tij  at 
a  temperature  T.  If  a  diffusive  flux  jn  arises  from 
inhomogeneous  doping,  in  equilibrium  it  must 


be  balanced  by  the  opposh-  Tut  flux  ;Dr 

;d  =  /dr  {••■ 

/dr  =  Hi  X  E;  jD  =  -  c,  ,1  c. J  (Hi) .  (45) 

Then,  employing  tlw  Einst.  in  relationshii 

Di  =  \iixkaT/ei  (46) 

we  obtain  for  the  electric  field  strength  E 

n,xE  =  (kff/ej)  x  grad  («,).  (47) 


In  eq  44-47,  D„  p,  and  e,-  are  the  diffusion 
coefficient,  the  mobility  and  the  electric  charge  of 


Figure  23  Static  electric  potential  at  the  surface  of 
freshwater  columnar  ice.  An  area  of  2.5x2. 5  cm2  is 
shown.  The  columns  lie  along  the  x-axis.  The  potential 
is  given  in  volts ;  T  =  -20  °C  (from  Petrenko  '1992). 


the  (th  charge  carrier  respectively;  k$  is  Boltz- 
man's  constant.  For  the  potential  difference  A<p 
between  two  points 


Atp  =  (cpi  -  (p2)  =  (A:b T/Ci)  x  ln(i(2l/»|1)  (48) 

where  n  1(-  and  Mj,-  are  the  concentration  at  points  1 
and  2.  Equation  48  gives  A(p  =  97  mV  for  Bjerrum 
defects  (L-  and  D-defects)  at  T~  263  K  and  «2// « 1/ 
=  10.  A  typical  magnitude  of  the  electric  field 
associate  I  with  Bjerrum  defects  in  ice  is  about  1- 
5  V/cm  (see  Fig.  23).  Since  A<p  depends  upon  a 
ratio  of  impurity  concentration,  such  an  electric 
field  exists  even  in  very  pure  single  crystals  of  ice. 

I  (Petrenko  1992)  have  studied  electric  poten¬ 
tial  on  tire  surface  of  natural  and  artificial  ice 
using  a  scanning  microprobe  connected  to  an 
electrometer  (Fig.  24).  The  central  element  of  the 
experimental  f-  r  is  an  X-Y  pen  recorder,  in 


which  a  pen  is  replaced  by  a  thin  stainless  steel 
electrode  with  a  ball  of  0.3  mm  diameter  at  its  end. 
A  computer  controls  the  position  of  the  electrode 
scanning  over  a  given  area  of  an  ice  sample.  Scan¬ 
ning  was  from  lxl  mm  to  5x5  cm,  with  resolution 
better  than  or  equal  to  ±0.5  mm.  The  pressure  on 
the  electrode  was  adjusted  so  finely  that  the  elec¬ 
trode  did  not  produce  visible  tracks  on  or  damage 
to  the  ice  surface.  The  velocity  of  the  mobile  elec¬ 
trode  and  the  measuring  time  can  be  changed 
within  a  wide  range,  so  that  tire  electrode  is  able  to 
measure  a  correct  value  of  an  electrochemical  sur¬ 
face  potential.  I  found  that  an  electrode  velocity  of 
less  than  10  mm  s-1  and  a  measuring  time  longer 
than  10  ms  were  satisfactory.  The  potential  was 
amplified  by  an  electrometer  preamplifier  with 
input  impedance  Rjnput  - 1 0 * 2  li  and  then  stored  in 
the  computer.  The  smooth  surface  of  the  ice  sam¬ 
ples  was  prepared  using  a  microtome. 

An  example  of  tire  electric  potential  over  tire 
surface  of  freshwater  columnar  ice  is  shown  in 
Figure  23.  lire  ice  columns  lie  along  the  x-axis.  A 
comparison  of  the  electric  potential  pattern  with 
an  optical  image  of  the  ice  sample  placed  between 
two  crossed  polarizers  reveals  tlreir  similarity  or 
even  identity.  In  Figure  23  "rough"  and  "smooth" 
bands  belong  to  different  individual  columns.  The 
most  significant  gradient  of  potential,  and  hence 
electric  field  strength  E  -  grad(  <|>),  is  located  at  or 
near  grain  bounds/ es.  Generally,  the  distribution 
of  the  electric  potential  in  ice  follows  the  structural 
pattern  of  the  ice;  arrangement  of  grains,  columns 
and  boundaries  between  the  regions  of  doped  and 
normal  ice. 

Figure  25  depicts  the  electrical  potential  over 
the  surface  of  a  two-part  ice  sample.  One  part  was 
grown  from  very  pure  deionized  water  and  one 
was  grown  from  water  doped  with  101K  molecules 
enr3  of  KOH.  Those  two  halves,  with  ideal  flat 
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Figure  25.  Electric  potential  ok  the  surface  of  ice  sample 
assembled  of  pure  and  KOH-doped  ice.  The  steep  step  of 
the  potential  is  located  on  the  boundary  between  the  pure 
and  doped  ice.  A  "ditch  "  seen  in  the  potential  of  the  doped 
(more  positive)  half  corresponds  to  locations  of  grain 
boundaries.  Dimension  of  the  area  shown  is  1x1  cm;  T  = 
-19.5°C  (from  Petrenko  1992). 

and  smooth  surfaces,  were  pressed  together  for 
10  minutes  at  T  =  -5°C,  so  that  they  built  a  mono¬ 
lithic  block.  Doping  ice  with  KOH  provides  nega¬ 
tively  charged,  extrinsic  mobile  defects:  L-defects 
and  OH" ions  (Howe  and  Whitworth  1989).  Their 
diffusion  into  pure  ice  makes  tire  pure  part  of  the 
sample  negatively  charged  and  the  doped  part 
positively  charged.  That  is  precisely  what  we  can 
see  in  Figure  25.  At  temperatures  higher  than 
approximately  -15*C  for  pure  ice  and  -22°C  for 
saline  ice,  the  existence  of  a  conductive,  quasi¬ 
liquid  layer  on  the  ice  surface  masks  (or  screens) 
the  electric  potential  measured  with  this  tech¬ 
nique 

PHOTO-EMISSION  OF  ELECTRONS 
FROM  THE  ICE  SURFACE 

Nason  and  Fletcher  (1975)  studied  photo¬ 
emission  of  electrons  (extrinsic  photo-effect)  from 
the  surface  of  ice  and  water.  They  found  that  the 
threshold  of  photo-emission  for  both  ice  near 
-10°C  and  water  near  10°C  is  about  195±5  run. 
This  wavelength  corresponds  to  the  energy  of  6.3 
±0.2  eV.  Hie  emission  efficiency  for  liquid  water 
was  about  half  that  for  ice.  Photo-emission  from 
ice  decreases  steadily  with  decreasing  tempera¬ 
ture,  apparently  without  any  shift  in  the  thresh¬ 
old  value,  and  vanishes  below  approximately 
-160°C.  This  decrease  of  photo-emission  with 
temperature  was  not  related  to  the  decrease  of 
conductivity  of  ice,  since  the  same  temperature 
dependence  of  photo-emission  was  found  for  ice 
heavily  doped  with  ITF  (5%).  On  the  basis  of 


similarity  of  photo-emissions  from  ice  and  water, 
the  authors  concluded  that  a  quasi-liquid  layer 
exists  at  the  ice  surface  at  temperatures  T  > 
-160°C.  Since  the  threshold  is  extremely  sensitive 
to  the  value  of  the  surface  potential,  its  weak 
sensitivity  to  temperature  variations  implies,  as 
the  authors  believe,  that  the  electric  structure  (A^, 
the  screening  length)  of  the  surface  layer  of  ice 
and  water  remains  unchanged  within  the  whole 
temperature  range  (-160  to  10°C).  Nason  and 
Fletcher  believe  that  the  most  probable  source  of 
photo-electrons  is  the  process  of  ionization  of 
OH~  ions  in  the  surface  layer  of  ice  and  water. 
This  hypothesis,  obviously,  contradicts  the  fact 
that  a  change  in  pH  of  ice  from  1  to  14  does  not 
affect  photo-emission  significantly. 

At  the  present  time,  we  can  suggest  an  inter¬ 
pretation  that  is  different.  It  is  known  that  illumi¬ 
nating  ice  with  light  having  a  wavelength  X  £  190 
nm  excites  a  bulk  photo-conductivity  according 
the  following  reaction 

H20  +  hv  -4  H20*  -4  HzO+  +  e  (49) 

where  H20*  denotes  an  excited  water  molecule 
(Khusnatdinov  and  Petrenko  1992).  The  electrons 
produced  in  this  reaction  can  tunnel  out  of  the 
ice,  provided  that  its  surface  has  a  negative 
potential.  Since  tire  reaction  (eq  49)  begins  with 
the  excitation  of  a  water  molecule,  it  is  not  sur¬ 
prising  that  it  is  insensitive  to  the  change  in  pH. 
The  temperature  dependence  of  photo-emission 
(as  well  as  the  difference  in  its  values  for  ice  and 
water)  might  result  from  the  change  in  the  sur¬ 
face  electric  field,  since  this  field  causes  tunneling 
of  electrons  out  of  the  subsurface  layer  of  tire  ice. 

SURFACE  OPTICAL  ABSORPTION 
IN  INFRARED  REGIONS 

Interesting  results  about  die  molecular  struc¬ 
ture  of  the  ice  surface  were  obtained  most  recent¬ 
ly  from  the  study  of  infrared  optical  absorption 
spectra  (Bush  and  Devlin  1991,  Rowland  and 
Devlin  1991,  Callen  et  al.  1992,  Devlin  1992,  Nix¬ 
on  et  al.  1992,  Rowland  et  al.  1993).  Those  studies 
and  the  resulting  implications  were  supported  by 
computer  simulations  using  molecular  dynamics 
simulation  methods  (Bush  and  Devlin  1991). 

The  above-mentioned  measurements  are 
based  on  the  possibility  of  distinguishing  be¬ 
tween  lines  in  the  IR  spectra  corresponding  to  the 
stretching  and  bending  mode  oscillations  of  non- 
saturated  hydrogen  bonds  (i.e.,  dangling  bonds) 
at  the  ice  surface.  Tie  frequency  of  these  oscilla- 


tions  appears  to  be  extremely  sensitive  to  the  way 
in  which  a  given  water  molecule  is  bonded  to 
other  molecules — by  two  or  by  three  hydrogen 
bonds  (2-coordinated  or  3-coordinated  water 
molecules).  Molecules  in  the  bulk  are  4-coordin- 
ated.  Therefore,  a  comparison  between  the  inten¬ 
sities  of  2-coordinated  and  3-coordinated  mole¬ 
cules  allows  us  to  estimate  their  relative  concen¬ 
trations.  The  IR  absorption  spectra  observed  in 
the  experiments  are  better  interpreted  because  of 
the  high  sensitivity  of  the  eigen  frequencies  of 
dangling  bond  oscillations  to  the  isotopic  substi¬ 
tution  of  H  with  D.  Moreover,  these  frequencies 
are  quite  sensitive  to  the  process  of  adsorption  of 
various  gases.  Thus,  such  measurements  allow  us 
to  study  ion  adsorption  by  the  ice  surface  at  low 
temperatures. 

Specimens  of  ice  with  large  surface-to-volume 
ratios  were  produced  either  in  the  form  of  micro¬ 
scopic  ice  clusters  suspended  in  gas  or  by  vapor 
condensation  onto  a  substrate,  forming  thin  ice 
films.  In  experiments  by  Rowland  et  al.  (1993),  for 
example,  crystalline  ice  clusters  -25  nm  in  diame¬ 
ter  formed  when  mixtures  of  an  inert  carrier  gas 
containing  a  small  amount  of  water  vapor,  with  a 
gas/H20  ratio  of  -200/1,  were  introduced  rapid¬ 
ly  into  a  cooled  cluster  cell  to  pressures  of  33  kPa 
(V 3  bar)  at  temperatures  below  -120  K.  Increas¬ 
ing  temperature  and  load  pressures  create  larger 
clusters  with  consequently  smaller  surface-to- 
volume  ratios  and  reduced  relative  surface-mode 
band  strengths  (Rowland  et  al.  1991). 

Computer  simulations  using  molecular  dy¬ 
namics  methods  have  demonstrated  that  when 
water  molecules  collect  to  produce  a  small  clus¬ 
ter,  tire  following  phenomena  occur;  1)  there  are 
no  molecules  possessing  two  dangling  OH 
bonds;  2)  both  2-  and  3-coordinated  molecules 
form  on  the  surface  in  approximately  equal 
amounts,  provided  that  the  temperature  is  quite 
low  (-20  K);  and  3)  in  the  annealed  clusters,  the  3- 
coordinated  species  are  dominant.  The  concen¬ 
tration  of  dangling  bonds  is  small  («  1)  even  at 
these  low  temperatures.  Tne  results  of  spectral 
studies  performed  in  the  above  works  are  in 
agreement  with  the  results  of  computer  simula¬ 
tions. 

REGELATION  (REFREE2ING), 

SINTERING  AND  ADHESION 

In  retrospect,  the  studies  of  the  quasi-liquid 
layer  at  the  ice  surface  began  with  studies  of  ther¬ 
mo-mechanical  effects:  regelation,  sintering  and 


adhesion.  The  d  tta  on  the  physical  parameters  of 
the  ice  surface  layer,  obtained  during  those  stud¬ 
ies,  are  ambiguous  and  much  less  accurate  than 
those  obtained  by  modern  methods.  But  we 
should  keep  in  mind  that,  first,  it  was  these  phe¬ 
nomena  that  aroused  such  a  great  interest  in  the 
ice  surface.  Second,  these  phenomena  are  still  of 
great  practical  importance.  Therefore,  a  brief  ret¬ 
rospective  view  of  these  studies  is  presented  be¬ 
low. 

Faraday  (1650)  observed  an  interesting  phe¬ 
nomenon:  if  two  pieces  of  ice  are  brought  into 
slight  contact  at  0°C  or  a  higher  temperature,  they 
freeze  together.  This  discovery  initiated  a  long 
and  interesting  argument  between  Faraday  and 
Tyndall  on  one  side  and  ].  Thomson,  his  brother 
YV.  Thomson  (Lord  Kelvin)  and  Helmholtz  on  the 
other  side.  The  Thomsons  and  Helmholtz  held  to 
the  opinion  that  regelation  is  caused  by  melting 
under  pressure  in  the  contact  zone  between  two 
pieces  of  ice.  When  the  pressure  is  removed,  wa¬ 
ter  freezes  again,  bonding  the  pieces  together. 
Faraday  performed  several  simple  but  ingenious 
experiments  that  proved  the  erroneousness  of  the 
Thomsons'  hypothesis,  suggesting  a  regelation 
mechanism  that  is  very  close  to  the  modem  con¬ 
cept  of  this  process.  He  assumed  that  at  tempera¬ 
tures  close  to  the  melting  point,  a  thin  liquid  film 
always  exists  at  the  ice  surface.  When  two  blocks 
of  ice  touch,  this  film  first  forms  a  "neck"  connect¬ 
ing  the  two  blocks.  Then,  since  the  inner  part  of 
the  neck  gradually  turns  into  a  bulk  region  of  ice, 
it  freezes,  forming  a  solid  joint  between  the  two. 

Nakaya  and  Matsumoto  (1953)  did  some  ex¬ 
perimental  work  to  verify  the  existence  of  the 
liquid-like  layer  at  the  ice  surface.  They  suspend¬ 
ed  small  spheres  of  ice  (from  1.5  to  4  mm  in 
diameter),  prepared  from  distilled  water,  on  cot¬ 
ton  threads  in  a  pendulum-like  fashion.  One  of 
these  pendulums  was  moved  laterally,  bringing 
the  spheres  in  contact  with  a  very  small  force, 
estimated  by  the  authors  to  be  less  than  0.01  dyn, 
Then,  the  suspension  point  of  one  of  the  spheres 
was  displaced  laterally  until  the  spheres  broke 
free;  the  force  can  be  calculated  from  the  value  of 
the  angle  0  (with  respect  to  the  vertical)  at  which 
they  separated.  Nakaya  and  Matsumoto  also  test¬ 
ed  ice  spheres  containing  0.1%  NaCl.  hr  this  case 
we  would  expect  a  thicker  film  on  the  ice  surface. 
To  explain  their  results,  the  authors  had  to  as¬ 
sume  that  a  transitional  liquid-like  layer  exists  at 
tire  ice  surface  at  temperatures  below  0°C. 

Hosier  et  al.  (1957)  did  similar  but  more  accu¬ 
rately  controlled  experiments  in  a  temperature 


% 


:v-  - W  i.  - 


■;  >J>- .'.T -k,*. 


in' h'i". "  '■  I 


22 


-BO  -25 


-15  -10 

Temperature,  “C 


Figure  26.  Separation  force  as  a  function  of  tempera¬ 
ture:  A — water-vapor  saturated  atmosphere ;  B — dry 
atmosphere  (after  Hosier  et  al.  1957). 


range  from  -80  to  0°C  in  both  dry  and  water- 
saturated  atmospheres.  Between  -80  and  -40°C, 
they  observed  no  adhesion  between  spheres  (Fig. 
26).  Tlie  lowest  temperature  at  which  adhesion 
was  measurable  in  a  saturated  atmosphere  was 
about  -25°C,  whereas  in  a  dry  atmosphere  the  ad¬ 
hesion  practically  ceased  at  -3°C.  They  explained 
this  difference  in  terms  of  a  rapid  evaporation  of 
the  liquid-like  layer  from  the  ice  surface  in  a  dry 
air,  which  makes  it  very  thin.  Hosier  et  al.  also 
found  that  the  temperature  dependence  of  the 
force  at  which  separation  occurs  can  be  expressed 
in  terms  of  a  Boltzmann  equation  with  the  activa¬ 
tion  energy  of  0.46  eV. 

Jellinek  (1957a, b;  1959;  1960;  1962),  investigat¬ 
ing  the  adhesive  properties  of  ice,  conducted 
shear  and  tensile  experiments  with  ice  frozen  on¬ 
to  stainless  steel  and  optical  flat  quartz  plates  and 
various  polymers,  hi  his  tensile  experiments, 
only  cohesive  breaks  through  ice  were  observed. 
Shear  tests  resulted  in  adhesive  breaks  along  the 
interface,  the  strength  of  which  depended  on  the 
loading  rate,  surface  polishing,  temperature  and 
the  type  of  the  material  used.  A  temperature  de¬ 
pendence  of  the  shear  strength  for  the  case  of  ice- 
snow  confined  between  stainless  steel  disks  is 
shown  in  Figure  27. 

If  we  assume  that  the  shear  strength  in  shear- 
type  experiments  at  temperatures  above  -13°C  is 
determined  by  viscous  flow  of  the  quasi-liquid 
layer,  then 
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Figure  27.  Strength  as  a  function  of  temperature  for 
snow-ice  sandwiched  between  stainless  steel  disks,  ob¬ 
tained  by  shear.  Cross-sectional  area  1.54  cm2;  height  0. 2 
to  0.4  cm.  Adhesive  breaks  at  T  >  -13’C  and  cohesive 
breaks  af  T  <-23  ‘C.  Each  point  represents  an  average  of  at 
least  12  tests  (after  Jellinek  1959). 

43m  a  9.  (50) 

dv  L 

where  q  and  L  are  film  viscosity  and  thickness 
respectively.  The  ratio  of  r\/L  found  by  Jellinek  at 
T  =  -4.5°C  was  6.99X104  N  s/m  for  ice-stainless 
steel  and  1.5  xlO4  Ns/m  for  ice-quartz.  In  eq  50 
we  have  two  unknowns,  q  and  L.  If  the  value  of  L  = 
2.7x  10'9  m,  found  by  Barer  et  al.  (1977)  at  T  =  -4°C 
(see  Table  1),  is  used,  then  the  resulting  viscosity 
of  the  quasi-liquid  layer  between  ice  and  quartz 
will  be  q  =  0.41  Pa  s,  as  compared  to  0.017  Pa  s 
found  by  Barer  et  al.  at  even  lower  temperatures. 

Raraty  and  Taber  (1958)  did  some  adhesion  ex¬ 
periments.  Their  results  are  essentially  in  agree¬ 
ment  with  those  of  Jellinek. 

A  large  number  of  experimental  and  theoretical 
works  dealt  with  the  issue  of  sintering  of  ice.  Sin¬ 
tering  plays  a  major  role  in  the  process  of  snow 
densification,  which  transforms  falling  snow  into 
ice  in  mountain  glaciers  and  tire  polar  caps  of 
Greenland  and  Antarctica.  Several  physical  phe¬ 
nomena,  such  as  surface  and  bulk  diffusion,  sur¬ 
face  tension,  plastic  deformation  and  water  trans¬ 
fer  through  water  vapor,  contribute  to  the  sinter¬ 
ing.  The  relative  effect  of  each  of  the  mechanisms 
above  depends  largely  on  conditions  during  sin¬ 
tering:  temperature,  humidity,  impurity  composi¬ 
tion  and  size  of  ice  particles.  Because  of  this,  inter¬ 
preting  the  data  on  ice  sintering  is  usually  very 
difficult.  The  phenomena  of  small  particle  migra- 
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tion  (Itagaki  1967a)  on  the  ice  surface  and  regela¬ 
tion  are  also  similar  to  sintering.  Studies  of  these 
things,  although  contributing  significantly  to  the 
development  of  practical  applications,  were  of  lit¬ 
tle  use  for  the  determination  of  the  fundamental 
parameters  of  the  ice  surface  layer  owing  to  diffi¬ 
culties  in  interpretation.  As  examples  of  excellent 
studies  on  ice  sintering,  I  can  recommend  to  inter¬ 
ested  readers  the  following  papers:  Kingery 
(1960),  Hobbs  and  Mason  (19641),  Itagaki  (1967c) 
and  Maeno  and  Ebinuma  (1983). 


ADSORPTION  OF  GASES 
ON  THE  ICE  SURFACE 


When  gas  molecules  of  a  mass  m  strike  a  solid 
surface,  they  can  either  bounce  or  get  stuck  to  the 
surface  (be  entrapped)  for  a  certain  time.  The  latter 
case  is  called  an  adsorption  of  the  molecule,  If  the 
gas  molecules  are  in  thermal  equilibrium  with  the 
solid  surface,  the  equilibrium  surface  concen¬ 
tration  Cs  of  adsorbed  molecules  will  depend  on 
tire  flux  of  oncoming  molecules,  i.e.,  on  tire  gas 
pressure  P  and  temperature  T  and  on  the  average 
time  that  a  molecule  spends  on  the  surface. 

Langmuir  (1916)  assumed  that  the  maximum 
concentration  of  the  adsorbed  molecules  Cm  corre¬ 
sponds  merely  to  a  mono-molecular  layer.  There¬ 
fore,  a  molecule  coming  from  a  gas  and  colliding 
with  an  already  adsorbed  molecule  must  instantly 
bounce  off  tire  surface.  On  the  basis  of  this  as¬ 
sumption,  simple  kinetic  considerations  lead  to 
the  so-called  Langmuir's  adsorption  isotherm 


Vk  -  kiP 
1/m  'l+kiP 


(51) 


where  l2,  is  the  actual  volume  of  the  adsorbed  gas 
(which  corresponds  to  Cs)  and  Vm  (which  corre¬ 
sponds  to  Cm)  is  the  volume  of  tire  adsorbed  gas 
when  the  entire  adsorbent  surface  is  covered  with 
a  monolayer. 

ki  - - I - (52) 

(Intti  hiT)1/7  Cm 


Brunauer  et  al.  (1938)  extended  tire  result  of 
Langmuir  to  the  case  of  multilayer  adsorption. 
This  equation  is  known  as  a  BET  equation.  They 
found  that 


P _ _  1  + 

V^P-Po)  VMk2  VMk2i\> 


(53) 


where 


k2  ~  exp 


(eo-e'I 
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Ea  and  Ea  are  the  heats  of  adsorption  (i.e.,  the 
amounts  of  heat  liberated  when  a  molecule  is  ad¬ 
sorbed)  corresponding  to  the  first  and  all  succeed¬ 
ing  layers;  P0  is  the  saturation  pressure  of  a  gas 
over  a  plane  surface  at  temperature  T. 

As  we  can  see  from  eq  53,  a  plot  of  P/ Vs  ( P  -  PQ) 
versus  P/Pa  must  be  a  straight  line  with  a  slope  of 
(k2  -  1)/ V'm  k2  and  an  intercept  with  a  horizontal 
axis  of  1  /  V'm  In  this  way  constants  VM  and  k2 
can  be  determined;  then,  the  surface  area  of  the 
solid  can  be  calculated.  We  can  also  find  Ea,  pro¬ 
vided  that  the  value  of  Ea  is  known  (it  is  usually 
taken  to  be  equal  to  the  latent  heat  of  evaporation 
of  the  adsorbate  in  a  liquid  form).  Readers  interest¬ 
ed  in  tire  theories  of  adsorption  can  find  them  in 
the  book  by  Steele  (1974),  for  example. 

Experimental  studies  of  gas  adsorption  on  the 
ice  surface  are  of  a  great  practical  importance  for 
the  physics  of  the  atmosphere,  meteorology,  cli¬ 
matology  and  astrophysics.  In  addition,  such  stud¬ 
ies  contribute  to  the  understanding  of  tire  ice  sur¬ 
face  structure,  since  adsorption  is  quite  sensitive  to 
both  surface  structure  and  tire  presence  of  various 
defects  at  tire  surface,  as  well  as  to  the  local  fields 
(from  tire  mohcular  dipoles,  for  instance)  existing 
near  the  surface. 

Gas  adsorption  by  ice  has  beexr  studied  in  a 
wide  temperature  range:  from  77  to  273  K.  The  ad¬ 
sorption  of  gaseous  nitrogen  was  studied  by  Ad¬ 
amson  et  al.  (1967),  Ocampo  and  Klinger  (1983a) 
and  Schmitt  et  al.  (1987);  the  adsorption  of  CO2 
was  studied  by  Ocampo  and  Klinger  (1983);  tire 
adsorption  of  hydrogen  was  examined  by  Smolu- 
clrowski  (1983;  theoretical  study);  and  Orem  and 
Adamson  (1969)  studied  the  adsorption  of  non¬ 
polar  hydrocarbon  vapors  (n-alkanes).  In  these  ex¬ 
periments,  the  surfaces  of  amorphous,  cubical  and 
hexagonal  ice,  cleaved  at  77  K.  and  annealed  at 
higher  femperatures,  and  natural  snow  were  in¬ 
vestigated. 

Cubic  ice,  produced  by  the  direct  condensation 
of  water  vapor  oirto  the  walls  of  a  vessel  at  a  liquid 
nitrogen  temperature,  Iras  the  maximum  net  sur¬ 
face  (11.8  nr2/g  [Adanrson  et  al,  1967]).  This  sur¬ 
face  was  reduced  after  the  ice  had  been  annealed, 
and  was  about  1  nr2/ g  in  the  case  of  natural  snow. 
The  value  of  the  adsorption  heat  found  exceeds  the 
heat  of  evaporation:  2700±20  J  mol-1  at  77  K  for  the 
adsorption  of  N2  reported  by  Schmitt  et  al.  (1987); 
a  similar  value  of  2512±250  J/mol  was  reported  by 
Adamson  et  al.  (1967).  The  adsorption  heat  is  gen- 
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erally  reduced  after  annealing  (down  to  1067  J 
mol-1  after  4  days  [Ocampo  and  Klinger  1983a]). 

The  results  of  studies  performed  by  Orem  and 
Adamson  (1969)  showed  that  there  exists  a  signif¬ 
icant  difference  in  the  character  of  adsorbent- 
adsorbate  interaction  above  and  below  -35°C.  Be¬ 
tween  -35  and  -196°C,  ice  acts  as  a  low-energy 
adsorbent  and  the  adsorption  energy  is  associated 
with  dispersion  interactions,  with  no  effect  of  lo¬ 
cal  dipole  fields.  However,  at  temperatures  above 
-35°C,  the  ice  surface  acquires  an  adsorption 
capability  close  to  that  of  a  water  surface.  These 
observations  conform  to  the  model  of  a  quasi¬ 
liquid  film  at  the  ice  surface.  A  similar  temp°ra- 
ture  transformation  in  tire  adsorption  process,  _  ,<t 
at  a  lower  temperature  (-40°C),  was  also  noted  by 
Ocampo  and  Klinger  (1983a).  They  also  four  J  a 
considerable  decrease  in  the  entropy  of  the  ad¬ 
sorbed  gas  at  T  >  -40°C.  Such  an  entropy  decrease 
could  be  explained  in  terms  of  the  ordering  of 
adsorbate  molecules  in  a  strong  electric  field  of 
ordered  water  molecules  at  the  ice  surface.  As  we 
have  seen  previously,  tins  also  conforms  to  the 
model  of  the  quasi-liquid  layer. 

SURFACE  ENERGY 

A  net  surface  free  energy  or  simply  surface 
energy  is  defined  as  the  work  required  to  increase 
the  interface  area  by  a  unit  area  (at  a  constant 
temperature).  We  shall  need,  for  a  theoretical 
analysis,  the  energies  of  the  interfaces  betw- 
solid  ice  (S),  liquid  phase  (L)  and  water  vapor  t  V)- 
We  denote  them  as  Ysl-  Yvs  and  Yvl  for  solid/ 
liquid,  vapor/ solid  and  vapor /liquid  interfaces 
respectively.  All  of  these  quantities  are  tempera¬ 
ture  dependent;  it  does  not  make  sense  to  gi\  e  a 
value  toy  without  indicating  the  temperature.  The 
best  known  quantity  is  the  most  easily  measui  i, 
Yvl-  which  can  be  determined  by  a  variety  of  <- 
perimental  techniques  based  on  the  measurem  .at 


of  the  surface  tension  of  water.  (These  techniques 
include  measurement  of  the  height  of  water  rise  in 
a  capillary  and  measurement  of  the  wavelength  of 
surface  waves,  excited  by  an  ultrasonic  generator, 
for  example.)  At  the  melting  temperature,  Yvl  * 
75.7  mj  rrr2  and  increases  with  increasing  temper¬ 
ature.  The  two  other  energies  Ysl  and  Yvs  cannot 
be  measured  quite  so  easily,  and  to  determine 
them  a  certain  ingenuity  is  required. 

We  shall  not  describe  here  the  early  and  less 
accurate  works,  and  shall  proceed  to  an  interest¬ 
ing  technique  suggested  by  Ketcham  and  Hobbs 
(1969).  That  technique  was  also  used  later  by  Har¬ 
dy  (1977).  The  technique  is  based  on  the  measure¬ 
ment  of  equilibrium  angles  at  tire  interfaces  be¬ 
tween  two  ice  seeds  and  water  vapor  or  water 
(Fig.  28).  As  can  be  seen  from  this  figure,  the  triple 
junction  in  a  grain-boundary  groove  is  subjected 
to  three  forces  of  surface  tension.  Therefore,  for 
equilibrium,  the  following  conditions  must  be  sat¬ 
isfied 


Ygb  =  2Yvs  cos(flVS/2) 

(55) 

Ygb  =  2Yls  cos(iYls/2)  . 

(56) 

Since  there  are  three  unknowns  Ygb  -  Yls  and  Yvs / 
we  need  a  third  equation  in  addition  to  eq  52  and 
53.  This  equation  can  be  the  condition  of  equilibri¬ 
um  among  vapor,  liquid  and  ice  in  the  case  of 
incomplete  wetting  (0  <r  (j>  <  n).  In  this  case 

Yvs  =  Yls  + Yvl  cost]).  (57) 

Using  this  technique  for  determination  of  the 
surface  energies,  Ketcham  and  Hobbs  (1969) 
found  that  at  0°C 

Yvs  -  (109  ±  3)  mj  nr2 

Yls  =  (33  ±  3)  mj  nr2  (58) 
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a.  Vapor/solid  interface. 


Liquid 

"M  6,s  t  + 

\\  // 


Ice  grain  1 


Ygb 


Ice  grain  2 


b.  Solid/liquid  interface. 


Figure  28.  Schematic  of  grain-boundary  grooves  ( after  Ketcham  and  Hobbs 
1969). 
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Ygb  =  (65  ±  3)  mj  irr2. 

The  reader  can  find  the  data  on  the  values  of 
the  angle  (J)  in  the  papers  by  Knight  (1966,  1971). 
Ketcham  arid  Hobbs  also  calculated  that  the 
anisotropy  of  Yvs  and  Yls  with  respect  to  orienta¬ 
tion  is  probably  small  and  falls  within  the  experi¬ 
mental  error.  The  value  of  ygl,  depends,  of  course, 
on  the  grain  boundary  type  and  must  tend  to  zero 
as  the  angle  of  misorientation  between  the  seeds 
tends  to  zero.  Hardy  (1977)  using  a  similar  tech¬ 
nique,  found  in  the  case  of  water/ice  interface  a 
more  exact  value 

Yls  =  (29.1  ±  0.8)  m]  nr2.  (59) 

Recently,  Makkonen  (in  press)  concluded  from 
results  of  measurements  of  <j>  (37°)  that  Yls  =  77  mj 
m-2  for  ice  at  -25°C. 

Finally,  Kuroda  and  Lacmann  (1982),  from  the 
number  of  broken  bonds  at  the  ice  surface  and  the 
value  of  heal  of  sublimation  of  8.5xl0_2tl  J  mole¬ 
cule-1,  estimated  Yvs  on  the  basal  (0001)  and  pris¬ 
matic  (1010)  faces  to  be 

Yvs  (0001)  -  121  mj  nr2 

Yvs  (1010)  =  128  mj  m"2  (60) 


SUMMARY  OF 

MAIN  EXPERIMENTAL  RESULTS 

Before  we  proceed  to  the  review  of  theoretical 
models  developed  to  explain  the  structure  and 
physical  properties  of  the  ice  surface,  let  us  sum¬ 
marize  the  major  experimental  results  so  that  we 
have  a  starting  point  for  comparing  them  with 
theories.  The  reader  undoubtedly  already  antici¬ 
pates  tire  difficulty  of  that  task.  The  range  of  dis¬ 
crepancies  in  the  results  of  different  authors  is  too 
large,  especially  when  different  experimental 
techniques  are  used. 

Nonetheless,  let  us  try  to  establish  the  facts 
that  either  conform  to  the  majority  of  results  or  to 
the  most  reliable  of  them.  I  shall  also  try  to  point 
out  major  discrepancies. 

All  results  can  be  classified  into  several  issues: 

1.  The  presence  or  absence  at  the  ice  surface  of 
a  special  layer,  which  has  physical  properties  dif¬ 
ferent  from  those  of  the  bulk  ice. 

2.  The  layer  thickness  and  its  temperature  do¬ 
main  of  existence. 

3.  The  molecular  structure  of  this  layer. 

4.  A  comparison  between  measured  properties 


of  the  ice  surface  layer  (diffusion,  electric  conduc¬ 
tivity,  refraction  coefficient,  etc.)  and  those  of 
water  and  ice. 

After  this  is  done,  we  shall  see  whether  it  is 
possible  to  assemble  some  consistent  picture  of 
the  ice  surface  layer. 

Existence  of  a  special  layer 

The  greatest  consistency  is  achieved  here.  Ex¬ 
perimental  papers  published  over  the  past  20 
years  never  address  the  question  of  the  existence 
of  the  special  layer.  All  authors  admit  to  its  exist¬ 
ence.  Let  us  join  them,  keeping  in  mind  an  im¬ 
pressive  listing  of  experimental  evidence  consid¬ 
ered  above. 

Thickness  of  the  layer  and 
the  domain  of  its  existence 

Since  in  all  cases,  a  strong  temperature  depen¬ 
dence  of  the  surface  properties  was  observed,  we 
cannot  discuss  the  question  of  the  layer  thickness 
without  referring  to  its  temperature  dependence. 

Figure  29  attempts  to  summarize  all  tire  mea¬ 
surements  of  the  thickness  of  the  surface  ice  layer. 
In  those  cases  when  a  temperature  dependence  of 
some  property  of  this  layer  (e.g.,  surface 
conductivity)  rather  than  the  absolute  thickness 
was  measured,  this  temperature  dependence, 
normalized  by  the  value  at  some  temperature  7'u, 
is  plotted. 

At  first  glance,  the  results  presented  in  Figure 
29  are  awesome,  so  great  are  the  discrepancies. 
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Figure  23.  Temperature  dependencies  of  quasi-liquid 
layer  thickness  taken  from  the  papers  cited  in  this 
review.  Capital  letters  at  the  curves  stand  for  first  letters 
of  the  authors'  names. 
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Thus,  the  layer  thicknesses  found  by  different 
authors  differ  by  one  to  two  orders  of  magnitude 
(at  -10°C,  L  «  10  A  [Elbaum  et  al.  1993]  and  L  « 
1000  A  [Golecki  and  Jaccar^  1978j,  for  instance). 

The  width  of  the  temperature  domain  of  exist¬ 
ence  also  differs  by  two  orders  of  magnitude  (AT 
«  1°C  [Elbaum  et  al.  1993]  and  AT  Z  100°C  [Mi- 
zuno  and  Hanafusa  1987]).  Taking  into  account 
the  thoroughness  of  the  measurements  and  con¬ 
trol  of  the  conditions  in  these  experiments,  I  feel 
that  it  is  most  unlikely  that  the  described  discrep¬ 
ancies  are  attributable  to  errors  in  measurement 
or  some  other  error.  On  the  contrary,  I  believe 
drat  all  measurements  were  correct,  but  absolute¬ 
ly  different  layers  on  the  surface  of  the  same  ice  Ih 
were  measured.  In  other  words,  the  surface  layer 
must  consist  of  several  layers,  each  possessing 
specific  physical  properties,  different  from  that  of 
the  bulk. 

In  fact,  optical  investigations,  based  on  the 
measurement  of  intensity  of  the  light  reflected 
from  the  surface,  are  sensitive  to  the  molecular 
density  of  the  layer.  Such  investigations  indicate 
the  presence  of  an  optically  dense  film  at  the  ice 
surface  in  the  temperature  range  from  0  to 
approximately  -A  to  -5°C.  The  refractive  index  of 
this  film  n  =  1.330,  determined  by  Furukawa  et  al. 
(1987a,b),  is  very  dose  to  the  refractive  index  of 
water  at  0°C  n  -  1.333. 

The  thickness  of  this  thinnest,  "external"  dense 
film  depends  not  only  on  the  temperature,  but 
also  on  the  type  of  crystallographic  face,  as  well 
on  the  presence  or  an  absence  of  air  (Elbaum  et  al. 
1993).  In  approximately  the  same  temperature 
interval  (-2.5  to  0°C),  Kouchi  et  al.  (1987)  ob¬ 
served  by  means  of  X-ray  diffraction  a  film  with 
no  long-range  order,  which  is  characteristic  of  a 
liquid  phase. 

Thus,  we  can  conclude  that  in  tne  mperature 
range  from  0°G  tr  "minur.  several  -b.c-es  O' 
sius,"  a  film  similui  in  structure  ana  density  to 
water  exists  on  the  surface  of  pure  ice. 

A  layer  with  increased  conductivity  can  be  de¬ 
tected  at  the  ice  surface  at  much  lower  tempera¬ 
tures  (-25  to  -35°C).  Its  maximum  thickness  is 
limited  by  the  value  of  the  screening  length.  An 
idea  of  the  loiter  can  be  obtained  from  the  work 
by  Chrzanowski  (1988):  L  =  0.5  |i.m  at  T  =  -53°C. 
Charge  carrier  concentration  and  conductivity  in¬ 
crease  within  this  layer  a-'-  they  approach  tne  sur¬ 
face.  Estimates  made  in  the  section  on  the  surface 
charge  ar  d  Lie  surface  potential  indicate  that  the 
concentration  of  ions  and  Bjerrum  defects  nea, 
the  surface  can  reach  a  value  comparable  to  the 


molecular  concentration  (see  eq  14  and  15).  Such 
a  large  concentration  of  defects  itself  can  tie  suffi¬ 
cient  reason  for  a  disordering  of  the  structure. 

A  layer  with  an  abnormally  large  amplitude  of 
thermal  vibrations  is  observed  at  even  lower  tem¬ 
peratures  (=  -50°C  [Golecki  and  Jaccard  1978]). 
The  thickness  of  this  layer  at  T  «  -1  °C  reaches  the 
value  of  1000  A. 

Even  larger  is  the  temperature  domain  (T  2: 
-100°C)  of  a  layer  detected  by  means  of  NMR 
(Mizuno  and  Hanafusa  1987).  This  layer  can  be 
characterized  as  having  a  spin-lattice  interaction 
different  from  that  of  the  bulk,  or  as  a  layer  with 
increased  self-diffusion  coefficient  and  rate  of 
molecule  reorientation.  The  thickness  of  this  lay¬ 
er  at  the  ice/  VerosU  interface  is  11 A  at  T  =  -12°C 
(Barer  etal.  1977). 

Molecular  swuchire  of  the  layer 

We  have  sever;  ,  sources  of  information  about 
Ore  surface  layer c  -ructure.  First,  as  was  shown  by 
Kouchi  et  al.  (1987),  in  the  temperature  interval 
0  to  -2.5°C,  there  is  no  long-range  order  at  the  ice 
surface,  which  is  characteristic  of  liquids.  The  dis¬ 
ordering  probably  remains  at  lower  tempera¬ 
tures,  but  the  technique  used  does  not  have  the 
r  acessary  sensitivity  to  detect  this.  An  interesting 
fact  is  the'  the  thickness  of  the  "liquid"  layer, 
which  is  measured  by  X-n.y  diffraction,  appears 
to  be  the  same  for  (0001)  and  (0110)  surfaces, 
which  contradicts  ellipsometrical  measurements. 
Therefore,  we  cannot  rule  out  the  possibility  that 
in  these  investigations  different  layers  were  ob¬ 
served. 

Golecki  and  Jaccard  (19/7,  1978),  who  used  a 
proton  channeling  technique,  observed  the  pres¬ 
ence  of  a  disordered  layer  down  to  -50°C.  ThL 
layer,  with  r  'Sickness  of  ■’l'out  1000  A  at  ‘ 
was  interpreted  as  havinu  .  greater  enq 
tude  of  the  thermal  vibrations  or  oxygen  aton. 
rather  than  having  an  amorpl'  our  rtructurc. 

Additional  indirect  evdence  of  tire  dir  jo. 
ing  in  tire  ice  surface  lay;  .  orncs,  as  vvv  !e  seen 
in  the  section  ;n  surface  conductivity,  from  a 
comparison  of  K.4R  data  with  the  data  on  Jre 
values  of  the  concentrations .  ind  the  mobilities  of 
charge  carriers  vsee  eq  11-11).  Wr>  found  that  the 
concentration  of  defective  hydrogen  bonds  (Ejer- 
rum  defects)  equals  2xl028  m"3  and  is  comparable 
to  tiie  molecular  concentration  of  3x1 023  m~3.  Ob¬ 
viously,  in  this  ease  it  is  impossible  to  preserve  an 
ordered  structure  interconnected  by  hydrogen 
bonds  corresponding  io  ice  I|v  Moreover,  since 
the  activation  energy  of  the  L-defect  mobility 
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(0.23  eV)  is  dose  to  that  of  <*=  G&  (0,29  eV,  Fig. 
9),  then  according  to  t_q  11  such  high  Bjerrum 
<?  ..  concentration  and  consequent  inevitable 
disci  dering  will  persist  down  to  very  low  tem¬ 
peratures  (*60°C).  We  should  keep  in  mind,  how¬ 
ever,  that  since  the  amplitude  of  the  NMR  signal 
decreases  with  decreasing  temperature,  the  thick¬ 
ness  'f  tills  la,  .«  must  also  decrease.  Finally,  let 
us  remember  that  tire  results  of  experiments  on 
the  surface  charge  density  and  gas  adsorption 
imply  that  a  fraction  of  water  molecules  in  the 
subsurface  layer  (at  least  within  a  single  mono- 
molecular  layer)  must  have  their  dipole  moments 
oriented  normally  to  the  surface.  Such  orientation 
is  also  an  ordering  of  some  kind. 


Physical  properties  of  the  surface  layer 
as  compared  with  water  and  ice 

As  we  discussed  earlier,  certain  physical  prop¬ 
erties  of  th“  surface  ice  layer  are  intermediate 
between  water  and  ice.  This  holds  for  viscosity 
and  self-diffusion  coefficient,  for  instance.  But 
there  are  two  remarkable  exceptions  to  this  rule, 
First,  the  optical  density  of  the  quasi-liquid  film 
practically  coincides  with  that  of  water  at  0°C. 
Second,  the  net  electric  conductivity  of  the  ice 
surface  exceeds  by  many  orders  of  magnitude  the 
electric  conductivities  of  both  ice  and  water  and 
corresponds  to  the  upper  boundary  of  conductiv¬ 
ity  of  superionic  materials. 
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PART  II.  THEORETICAL  MODELS  OF  THE  ICE  SURFACE 


We  have  seen  in  Part  I  how  complicated  the 
structure  of  the  ice  surface  is  and  what  diverse 
physical  properties  it  possesses.  It  is  clear,  there¬ 
fore,  how  difficult  it  is  to  create  a  theoretical  model 
describing  such  physical  objects  as  the  surface  lay¬ 
er  of  ice,  and  nobody  has  yet  been  able  to  devise 
one.  Theoretical  models  do  exist,  although  there 
are  very  few,  that  in  the  best  case  can  describe 
satisfactorily  only  one  of  the  phenomena  associat¬ 
ed  with  the  ice  surface  in  a  limited  temperature 
range.  However,  probably  all  the  physical  reasons 
that  make  the  ice  surface  as  it  is  were  analyzed  in 
these  works.  Therefore,  it  seems  reasonable,  in 
spite  of  a  very  limited  application  of  existing  theo¬ 
ries,  to  consider  and  discuss  them  briefly. 

The  most  simple  and  clear  is  the  phenomeno¬ 
logical  model  based  on  the  idea  of  a  possible  gain 
in  the  free  energy  of  tire  ice/vapor  interface  result¬ 
ing  from  the  formation  of  a  liquid  film  on  the  ice 
surface.  Obviously,  this  is  possible  only  if 


Tvs>  Yls  +  Yvl,  (61) 

Within  the  framework  of  this  theory,  the  sur¬ 
face  layer  of  ice  is  considered  to  be  simple,  ordi¬ 
nary  water.  If  we  denote  the  change  in  the  Gibbs 
free  energy  per  unit  area  resulting  from  melting  as 
AG,  obviously  in  thermodynamical  equilibrium 

(WjrpLAp(d)=0  (62) 

where  p  is  the  chemical  potential,  pL  is  the  water 
density  and  d  is  the  water  film  thickness.  The  wa¬ 
ter  film  thickness  d  enters  into  AG  twice 

AG  -  (pL  SL  AT)  d  +  Ay{d)  (63) 

where  AT  =  Tm  -  T,  Sl  is  the  entropy  of  fusion  per 
unit  mass  and  Ay  is  the  change  in  the  free  energy  of 
tiie  surfaces 


Ay(d)  =  Yvi  +  Ysl  -  Tvs-  (64) 

Ay  depends  on  the  film's  thickness  implicitly  be¬ 
cause,  at  thickness  values  comparable  to  the  inter¬ 
atomic  distances,  two  interfaces  begin  to  "sense" 
the  presence  of  each  other.  Usually,  it  is  assumed 
that 


A  fid)  -  A  y(°°)  [l  --  (~)2 


(65) 


where  r0  is  a  constant  on  the  order  of  interatomic 
distance  (see  for  example  Dash  [1991]  and  Takagi 
[1990]).  A  solution  of  the  system  of  eq  62-65  gives 
for  the  water  film  thickness 


2rp  A  y(qq)  \lfi 
,PlSl  (Tm  -  T)  j 


(66) 


The  reader  can  find  examples  of  these  calcula¬ 
tions  in  the  paper  by  Dash  (1991),  who  also  consid¬ 
ered  an  ice/solid  interface  and  took  into  account  a 
small  correction  due  to  thermomolecular  pressure 
acting  on  the  melted  layer;  also  in  the  paper  by 
Takagi  (1990),  who  incorporated  corrections  due 
to  surface  curvature.  Figure  30  shows  the  results  of 
such  a  calculation. 

Fukuta  (1987)  and  Makkonen  (in  press)  sug¬ 
gested  that  tine  surface  can  melt  because  of  the 
pressure  applied  to  the  subsurface  layer  by  the 
bulk.  The  nature  of  this  pressure  becomes  clear  if 
we  recall  that  a  molecule  on  the  surface  is  subject¬ 
ed  to  the  attraction  forces  from  the  rest  of  the 
molecules  and  the  net  force  is  directed  inside  the 
crystal.  Since  tire  pressure  lowers  the  melting  tem¬ 
perature  of  ice,  the  surface  melts.  Fukuta,  analyz¬ 
ing  the  intermolecular  forces  and  using  the  depen¬ 
dence  Tm(P),  valid  for  microscopic  quantities  of 
water  and  ice,  estimated  the  thickness  of  tire  liquid 
film  as 


d  -  1-5  nm 
VT^T  ' 


(67) 


Tm-T 

Figure  30.  Calculated  temperature  dependence  of  a  thick¬ 
ness  of  liquid  film  on  an  ice  surface  (after  Takagi  1990), 
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Makkonen  (in  press),  who  found  experimentally 
that  cq  61  is  not  valid,  predicted 


d  = 


3.5  nm  °C 


(68) 


(Tm-T) 

Within  the  context  of  Fukuta's  model,  the  liquid 
film  must  vanish  at  temperatures  T  <  -22 °C, 
where  the  pressure  just  transforms  ice  I},  into  its 
high  pressure  forms  (In,  Im  and  IIX)  rather  than 
into  liquid. 

Thus,  the  models  just  described  could  have 
been  applied  to  the  optical  observations  of  a  dense 
liquid  film  on  tire  ice  surface  close  to  the  melting 
point,  although  the  estimation  using  eq  67  gives  a 
much  thinner  thickness  compared  to  experimental 
data  (for  example  at  T  =  -1  °C;  see  Fig.  29).  hr  addi¬ 
tion,  these  phenomenological  models  can  give 
nothing  more  hr  terms  of  properties  of  the  liquid 
film  than  the  properties  of  ordinary  water.  And 
these  properties  do  not  coincide  at  all  with  the 
properties  of  tire  liquid-like  layer  of  ice,  as  we 
have  seen  previously  (electrical  conductivity  and 
self-diffusion  coefficient,  for  instance). 

During  a  rather  long  period,  a  model  proposed 
by  Fletcher  has  attracted  (and  still  continues  to 
attract)  considerable  attention  (see  Fletcher  1962, 
1963,  1968).  (Note  that  Fletcher  recently  returned 
to  the  problem  of  the  ice  surface  structure  at  much 
lower  temperatures  [see  Fletcher  1992].)  The  key 
idea  in  Fletcher's  (1968)  model  is  the  assumption 
that,  at  the  ice  surface,  a  significant  fraction  of 
molecules  (denoted  as  oj  is  oriented  with  their 
dipole  moments  (i.e.,  protons)  pointing  outwards, 
as  opposed  to  the  bulk  where  molecules  are  ori¬ 
ented  randomly.  This  would  result  in  a  large 
buildup  of  positive  polarization  charge  at  tire  sur¬ 
face.  Fletcher  came  to  his  conclusion  about  mole¬ 
cule  orientational  ordering  by  analyzing  the  inter¬ 
action  between  dipole  and  quadrupole  moments 
of  water  molecules  close  to  the  surface.  Fletcher 
estimated  the  energy  of  this  interaction  as  10_2°  J 
per  molecule,  which  exceeds  an  average  thermal 
energy  fcBT  <*  4xl0-21  J. 

As  we  have  seen  in  the  section  on  surface 
charge  and  surface  potential,  a  positive  electric 
charge  of  a  high  density,  which  is  unambiguously 
related  to  tire  water  molecule  dipole  ordering,  was 
indeed  discovered  on  the  ice  surface. 

All  other  predictions  of  Fletcher's  model  are 
consequences  of  the  molecule  ordering  in  the  sur¬ 
face  layer.  In  fact,  since  in  the  bulk  there  is  no 
ordering  (ac  =  x/i),  between  the  surface  and  the 
bulk,  there  should  be  a  transitional  layer  with  a 
large  concentration  of  defective  hydrogen  bonds 


allowing  molecules  to  reorient.  This  transitional 
disordered  layer,  according  to  Fletcher,  is  the  liq¬ 
uid-like  layer.  Fletcher  finds  the  equilibrium  thick¬ 
ness  of  this  layer  thermodynamically  by  minimiz¬ 
ing  the  Gibbs  free  energy  of  the  surface  and  incor¬ 
porating  the  following  components  of  this  energy: 
the  energy  associated  with  the  surface  orientation 
of  the  layer;  the  energy  associated  with  the  interac¬ 
tion  of  the  molecular  quadrupoles  with  the  gradi¬ 
ent  of  the  electric  field;  the  interaction  or  molecular 
dipoles  with  the  total  electric  field;  and  the  interfa¬ 
cial  energy  between  the  liquid-like  layer  and  tire 
ice.  In  addition,  a  "volume"  reduction  in  G  be¬ 
cause  of  the  melting  of  the  layer  and  a  lowering  of 
the  entropy  because  of  molecule  orientation  are 
apparently  also  considered.  An  unknown  fitting 
parameter  of  the  theory  is  £2,  the  energy  required 
to  form  an  L-defect  in  the  liquid-like  layer.  Finally, 
Fletcher  minimizes  the  free  energy  G  as  a  function 
of  ot0  and  the  layer  thickness  d.  I  deliberately  do 
not  present  the  mathematical  details  of  Fletcher's 
calculation  because  it  is  quite  cumbersome  and, 
more  importantly,  might  be  made  more  accurate 
in  the  future.  Nevertheless,  Fletcher's  analysis  of 
the  problem  and  tire  insight  that  he  has  achieved 
into  the  physical  nature  of  the  layer  are  absolutely 
remarkable.  Let  us  note  also  that  within  the  con¬ 
text  of  this  model,  an  unusually  high  electrical 
conductivity  of  tire  surface  layer  has  an  obvious 
explanation:  it  is  ascribable  to  a  high  concentration 
of  L-defects  in  the  liquid-like  layer.  Figure  31 


Figure  31.  Calculated  dependence  of  the  equilibrium 
thickness,  L,  of  the  liquid-like  layer  on  ice  for  various 
assumed  values  of  £2,  the  energy  of  formation  of  diffuse 
L-defects  in  the  liquid-like  layer  (after  Fletcher  1968). 
a— e2  =  6.25x1 0~3  eV;  b—e2  =  0.01 88  eV; 

c—e2  =  0.0625  eV;  d—e2  -  0.1875  eV. 
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Figure  32.  Schematic  representation  of  protonic  point  defects 
concentration,  n„  as  a  function  of  a  distance  from  the  ice 
surface ,  x;  d  is  the  thickness  of  a  superionic  layer.  The  nfniuo 
axis  is  logarithmic. 


0  <1 

shows  the  calculated  thickness  of  the  liquid-like 
layer  as  a  function  of  temperature,  Later,  in  his 
invited  review  paper,  Fletcher  (1973)  gave  an  ap¬ 
proximate  expression  for  the  thickness  d  of  the 
liquid-like  layer 

d( A)  =  (20  to  50)  -  25  logi0(273  -  T).  (69) 

Fletcher's  (1992)  recent  calculations  suggest 
that,  at  much  lower  temperatures,  tire  ice  surface  is 
reconstructed  by  an  orientation  of  surface  dipoles 
into  a  superlattice.  Such  a  reconstruction  should 
occur  on  the  (0001)  basal  plane  below  approxi¬ 
mately  30  K  and  on  the  (0110)  prismatic  plane 
below  approximately  70  K. 

On  the  basis  of  the  model  of  the  superionic  state 
of  ice  invented  by  Ryzhkin  (1985),  Ossipyan  and 
Petrenko  (1988)  suggested  a  purely  electrostatic 
reason  for  surface  melting  that  accounts  for  a  high 
surface  electric  conductivity  and  the  appearance  of 
a  rather  thick  (up  to  the  screening  length)  disor¬ 
dered  layer.  The  reasons  for  the  formation  of  the 
layer  are  as  follows.  As  was  shown  by  Ryzhkin 
(1985),  upon  the  increase  of  charge  carrier  concen¬ 
tration  above  a  certain  critical  value  (on  the  order 
of  ICH  to  10-4  of  the  molecular  concentration),  the 
activation  energy  of  ion  and  Bjerrum  defect  forma¬ 
tion  drastically  decreases,  leading  to  a  catastrophic 
growth  of  defect  concentration  up  to  0.1-1  of  the 
molecular  concentration.  This  is  equivalent  to  a 
phase  transition  into  a  superionic  state,  which  can 
not  have  a  regular  lattice  owing  to  a  high  defect 
concentration. 

If  the  surface  of  ice  is  charged  (which  is  the  case 
in  reality,  as  we  have  already  seen  in  the  section  on 
surface  charge  and  surface  potential),  this  charge 
will  be  screened  by  ions  and  Bjerrum  defects  com¬ 
ing  from  the  bulk,  as  shown  in  Figure  32.  By  taking 
into  account  the  transition  into  a  disordered 


superionic  state,  we  can  transform  a  smooth  curve 
of  the  screening  charge  concentration  distribution 
into  the  curve  with  a  tliin  superionic  layer  near  the 
surface.  Within  the  framework  of  this  model,  the 
high  surface  conductivity  of  ice,  the  formation  of 
the  quasi-liquid  disordered  layer  and  the  sharp 
interphase  boundary  (detectable  by  optical  meth¬ 
ods)  can  be  quantitatively  accounted  for.  To  my 
knowledge.  Dr.  Ryzhkin  is  currently  working  on 
the  improvement  of  this  model. 

There  exists  a  very  nontrivial  and  unusual  mod¬ 
el  of  the  structure  and  properties  of  the  ice  layer  at 
temperatures  close  to  the  melting  point.  Dzy- 
loshinskii  et  al.  (1961)  showed  that  a  fluid  that  is 
less  polarizable  than  its  solid  phase  is  attracted  to 
the  solid,  causing  a  liquid  film  to  grow  at  tine  sol¬ 
id/  vapor  interface.  Elbaum  and  Schick  (1991)  ap¬ 
plied  this  idea  to  the  ice  surface,  Figure  33  shows 
tine  result  of  their  calculations.  If  we  compare  this 


Figure  33.  Calculated  temperature  dependence  of  the 
liquid  film  on  ice  thickness  (after  Elbaum  and  Schick 
1991). 


figure  to  the  experimental  results  (Fig.  29),  it  be¬ 
comes  apparent  that  the  contribution  of  this 
mechanism  to  the  experimentally  measured 
thickness  of  the  quasi-liquid,  film  is  extremely 
small.  Kroes  (1992)  used  a  molecular  dynamics 


computer  simulation  to  study  the  surface  of  a 
small  cluster  of  water  molecules.  He  found  such 
ice  surface  characteristics  as  premelting  and  par¬ 
tial  polarization  of  molecular  dipoles. 
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